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SECTION  1 


INTRODUCTION 

This  report  deals  with  research  performed  at  UMIST  into  the 
electric  charge  that  is  transferred  when  vapor-grown  ice  crystals 
collide  and  separate  from  ice  particles.  The  objective  is  to 
understand  the  electrification  of  thunderstorms  in  which  it  is 
now  generally  accepted  that  the  electric  dipole  is  set-up  by 
charges  on  cloud  particles.  Our  laboratory  work  is  performed 
inside  a  cold-room,  in  which  a  water  droplet  and  ice  crystal 
cloud  may  be  grown  to  simulate  conditions  inside  a  thunderstorm. 
A  target  acts  like  a  graupel  pellet  (a  soft-hailstone)  which,  in 
a  thunderstorm,  falls  through  the  cloud  and  collides  with  super¬ 
cooled  water  droplets  and  ice  crystals.  This  ice-coated  target 
may  be  on  a  rotating  frame  which  sweeps  out  the  cloud,  or  it  may 
be  stationary,  in  which  case  the  cloud  is  drawn  past  the  target. 

Previous  work  in  UMIST  has  shown  that  at  cold  temperatures 
the  target,  simulating  a  graupel  pellet,  charges  negatively  with 
the  crystals  removing  the  equal  and  opposite  positive  charge.  In 
thunderstorms,  the  graupel  pellets  fall  and  lead  to  the  negative 
charge  center,  while  the  ice  crystals  are  carried  in  the  updraft 
to  form  the  upper  positive  charge  of  the  thunderstorm.  This 
upper-positive  lower-negative  configuration  has  been  well  estab¬ 
lished  as  being  typical  of  summer-type  thunderstorms  throughout 
the  world.  The  laboratory  work  has  also  shown  that  at  tempera¬ 
tures  warmer  than  typically  -20C,  the  graupel  pellets  charge 
positively  and  fall  to  form  a  lower  positive  charge  region,  which 
has  also  been  observed  in  the  base  of  thunderstorms.  In  this 
case,  the  crystals  are  carried  in  the  updraft  and  reinforce  the 
negative  charge  region  formed  from  the  falling  graupel  particles. 


The  laboratory  work  has  shown  that  the  temperature  at  which  the 
charge  transfer  changes  sign  is  a  function  of  the  liquid  water  in 
the  cloud. 

Prior  to  the  present  work,  our  laboratory  experiments  had 
shown  that  the  charge  transfer  was  a  sensitive  function  of  the 
interaction  speed  between  the  two  particles.  Also,  the  ice 
crystal  size  was  important  in  that  the  charge  transfer  depended 
on  size  to  the  fourth  power .  This  strong  dependence  was  used  in 
ext^  .  lations  of  the  laboratory  results  to  thunderstorms,  in  an 
attempt  to  account  for  the  total  charge  transfer  observed  in  the 
real  situation.  However,  the  laboratory  cloud  was  limited  in 
extent,  and  it  was  possible  only  to  grow  crystals  of  100  pm 
diameter;  whereas,  in  nature,  crystals  have  time  to  grow  several 
hundred  microns.  Another  conclusion  was  that  the  magnitude  of 
the  charge  transfer  was  affected  by  the  presence  of  liquid  water 
in  the  cloud.  With  a  cloud  that  consisted  of  ice  crystals  alone, 
the  charge  transfer  was  weak  and  inadequate  to  explain  thunder¬ 
storm  electrification;  whereas,  when  water  was  present,  substan¬ 
tial  charge  was  transferred. 

The  present  work  extends  the  previous  studies,  thanks  to  the 
provision  of  a  new  cold-chamber  with  greater  height.  Inside  this 
chamber,  crystals  of  several  hundred  microns  have  now  been  grown, 
and  these  have  been  used  in  a  series  of  experiments  which  have 
provided  some  surprising  new  results. 

In  addition  to  the  thunderstorm  experiments,  some  high  speed 
experiments  have  been  performed  to  determine  whether  the  charge/ 
velocity  results  obtained  previously  may  be  extrapolated  to  high 
velocities.  These  studies  will  have  some  relevance  to  the  impact 
of  cloud  particles  with  airplanes  which  may  lead  to  the  build-up 
of  static  charges  on  the  airplane  skin,  possibly  leading  to  a 


lightning  strike. 

This  report  is  divided  into  sections  dealing  with  various 
aspects  of  the  work.  The  basic  experimental  set-up  is  described 
in  Section  2,  together  with  details  of  some  of  the  experimental 
techniques  used.  Additional  information  is  located  in  the  appen¬ 
dices.  Section  3  deals  with  a  novel  method  of  measuring  the 
liquid  water  content  in  the  laboratory  cloud,  which  may  have 
application  to  measurements  in  the  natural  environment.  Section 
4  investigates  the  enhanced  collection  efficiency  of  a  rime- 
coated  target,  simulating  a  graupel  pellet,  for  super-cooled 
water  droplets.  The  non-uniform  rime  deposits  disturb  the  air¬ 
flow  around  the  collector  leading  to  enhanced  collection.  Sec¬ 
tion  5  covers  measurements  of  the  collection  efficiency  of  grau¬ 
pel  for  ice  crystals.  The  theory  of  this  topic  is  difficult  and 
very  few  laboratory  measurements  have  been  made.  The  results  are 
important  to  the  electrification  work,  as  electric  charge  is 
transferred  only  when  ice  crystals  separate  from  an  ice  pellet. 
Section  6  covers  the  charge  transfer  experiments  themselves. 
Section  7  discusses  the  results,  draws  conclusions  and  covers  the 
future  objectives  for  the  continuation  of  these  studies. 
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SECTION  2 


INSTRUMENTATION 

The  charge  transfer  experiments  were  performed  in  a  cloud 
chamber  situated  within  a  large  cold-room.  The  cold-room  has 
dimensions  of  3m  x  2m  x  4m  and  has  an  effective  operating  tem¬ 
perature  of  0°C  to  -378C.  The  refrigeration  equipment  has  a 
total  cooling  power  of  5.6  kW.  The  cloud  chamber  consists  of  a 
large  metal  box  divided  into  two  smaller  chambers,  as  shown  in 
Figure  1.  The  size  of  the  lower  chamber  is  0.94m  x  1.5m  x  2m  and 
the  upper  chamber  has  dimensions  0.9m  x  1.5m  x  lm.  The  two 
chambers,  which  are  connected  by  a  trapdoor  situated  in  the  floor 
of  the  upper  chamber,  remain  essentially  independent  of  one 
another.  The  two  chambers  have  a  number  of  observation  windows 
which  are  shielded  to  prevent  any  build-up  of  charges  affecting 
the  experimental  apparatus.  Figure  1  also  shows  the  positions  of 
the  six  thermocouples  which  were  used  to  monitor  the  chamber 
temperature.  The  cloud  chamber  was  earthed  to  prevent  the  app¬ 
aratus  from  being  affected  by  radio  frequency  noise.  The  two 
chambers  were  not  quite  air-tight  so  that  air  from  the  cold-room 
could  replace  air  that  was  removed  during  the  experimental  runs. 

The  cloud  was  formed  by  introducing  steam  from  a  boiler 
through  a  hole  in  the  floor  of  the  lower  chamber.  The  droplet 
cloud  in  the  upper  chamber  was  produced  by  an  ultrasonic  nebuli- 
ser  situated  on  a  shelf  outside  the  cold-room.  The  power  to  the 
kettle  could  be  controlled  from  outside  the  cold-room,  thus 
enabling  control  of  the  liquid  water  content  of  the  cloud  which 
was  measured  with  the  hot-wire  device  described  in  Section  3. 
The  accuracy  of  this  device  was  approximately  0.02  g  m~^  and  was 
insensitive  to  the  presence  of  ice  crystals.  As  noted  by  Mossop 
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(1984)  the  drop  size  distribution  could  be  modified  by  changing 
the  size  of  the  water  vapor  inlet  nozzle.  The  nozzle  diameter 
used  here  was  25  mm;  at  a  liquid  water  content  of  1  g  m  ^  the 
drop  size  distribution  had  a  maximum  size  of  33  pm  and  a  modal 
diameter  of  12  pm.  When  the  cloud  liquid  water  content  was  about 
2  g  m  ^  there  was  a  vertical  temperature  gradient  of  l^C  m  *  in 
the  lower  chamber. 

The  ice  crystals  were  produced  by  seeding  the  cloud  with  a 
wire  which  had  been  dipped  in  liquid  nitrogen.  The  wire  was 
introduced  quickly  into  the  cloud  chamber,  which  nucleated  the 
cloud  and,  as  the  environment  was  supersaturated  with  respect  to 
ice,  the  ice  crystals  grew  at  the  expense  of  the  water  droplets 
and  could  be  seen  as  'dust'  in  a  light  beam.  The  maximum  crystal 
size  which  could  be  produced  in  this  was  was  about  130  pm. 
Larger  ice  crystals,  >600  pm,  were  produced  using  a  slightly 
different  method.  The  droplet  cloud  was  formed  in  both  chambers; 
in  the  upper  by  the  nebuliser  and  in  the  lower  by  the  steam 
generator.  As  the  nebuliser  produced  droplets  by  ultrasonic 
vibration  and  not  boiling,  a  cloud  was  formed  in  the  upper  cham¬ 
ber  without  any  subsequent  rise  in  the  environment  temperature 
which  lead  to  the  upper  chamber  being  about  2°C  colder  than  the 
lower.  The  upper  cloud  was  seeded  using  the  method  described 
above,  and  a  suction  pump  connected  via  a  tube  to  the  top  of  the 
chamber  provided  an  updraft  which  countered  the  terminal  fall 
velocity  of  the  growing  crystals.  The  droplet  supply  was  main¬ 
tained  and  the  crystals  grew  to  about  150-200  pm  in  diameter. 
After  about  two  minutes,  the  suction  pump  wa.:  switched  off  and 
the  trapdoor  connecting  the  two  chambers  was  opened.  The  crys¬ 
tals  fell  through  into  the  supersaturated  environment  of  the 
lower  chamber  and  continued  to  grow.  About  30  seconds  later,  the 
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crystals,  now  approximately  300  |jm  in  diameter,  appeared  in  a 
light  beam  at  the  floor  of  the  lower  chamber.  The  crystals  were 
slowly  lost  to  the  walls  and  floor  of  the  chamber,  and  lasted 
only  for  a  further  two  minutes.  At  this  time  the  crystal  size 
reached  about  600  pm.  The  final  size  and  number  concentration  of 
the  ice  crystals  was  controlled  by  the  number  produced  in  the 
upper  chamber.  An  initial  high  concentration  led  to  a  high 
concentration  of  smaller  crystals.  Conversely,  a  low  number  led 
to  a  low  number  of  larger  crystals;  the  largest  diameter  ob¬ 
tained  was  over  700  pm. 

The  nature  of  the  cloud  particles  in  the  experiments  was 
determined  by  the  continuous  formvar  replicator  vhich  is  des¬ 
cribed  in  Appendix  I.  The  cloud  particles  are  drawn  past  a  film 
covered  with  2%  formvar  in  chloroform,  the  particles  stick,  and 
the  solution  dries  to  leave  a  plastic  replica  of  the  particle. 
The  replicator  had  collection  efficiencies  of  unity  for  droplets 
and  ice  crystals  greater  in  diameter  than  10  and  20  microns 
respectively.  The  film  was  marked  with  water  resistant  ink  to 
provide  time  marks  throughout  the  cloud  development.  The  replica¬ 
tor  could  not  collect  ice  crystals  greater  than  150  pm  in  dia¬ 
meter  at  the  flow  speeds  used.  As  the  crystals  fell  near  the 
replicator  orifice  they  experienced  a  horizontal  force  due  to  the 
air  being  drawn  into  the  orifice.  However,  due  to  the  inertia 
and  drag  of  the  large  crystals,  they  were  only  displaced  slightly 
towards  the  orifice  in  the  time  taken  to  fall  through  the  flow. 
Hallett  collected  1mm  diameter  ice  crystals  with  such  a  replica¬ 
tor,  but  at  aircraft  velocities.  These  speeds  are  not  practic¬ 
able  in  the  laboratory,  as  the  film  becomes  clogged  with  water 
droplets.  In  order  to  collect  the  large  crystals,  a  special 
method  was  developed  following  the  work  of  Schaefer  (1962),  van 
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den  Hage  (1969)  and  Saunders  and  Wahab  (1973).  The  cloud  parti¬ 
cles  were  collected  on  specially  prepared  microscope  slides, 
which  were  exposed  to  the  cloud  for  a  known  time.  The  slides 
were  prepared  in  the  laboratory  using  the  following  method.  A  4% 
solution  of  formvar  in  chloroform  was  made  in  the  usual  way.  The 
slides  were  then  dipped  into  the  formvar,  the  excess  solution  was 


drained 

Off 

and 

the 

slides  were 

placed  in  a  desiccator 

and 

al lowed 

to 

dry . 

The 

desiccator 

prevents  contamination  of 

the 

slides  by  water  droplets  condensing  on  the  formvar  by  evaporative 
cooling.  When  the  slides  were  required,  they  were  placed  in 
another  desiccator  in  the  cold-room  and  cooled  to  the  ambient 
temperature.  A  tube  of  diameter  5  cm,  which  is  connected  to  a 
vacuum  pump,  protrudes  into  the  cloud-chamber  at  the  same  at  the 
same  height  as  the  targets.  When  the  pump  is  switched  on,  the 
cold  dry  slide  is  inserted  into  the  tube  perpendicular  to  the 
flow  and  is  exposed  to  the  flow  for  five  seconds  and  then  removed 
from  the  tube.  The  slide  is  then  sprayed  with  a  fine  mist  of 
chloroform  from  an  atomiser  (kept  at  the  ambient  cold-roam  tem¬ 
perature)  to  activate  the  dry  formvar  and  then  replaced  in  the 
cold-room  desiccator  to  dry.  The  slide  could  be  removed  from  the 
cold-room  for  analysis  under  a  microscope  after  about  three 
minutes.  Crystals  >800  pm  have  been  successfully  replicated  with 
these  slides.  Previous  methods  have  required  that  the  formvar 
solution  is  kept  at  about  -5°C,  which  means  that  the  bulk  solu¬ 
tion  is  very  susceptible  to  corruption  by  water  vapour.  If  this 
occurs,  the  formvar  solution  is  useless,  as  hair-like  growth 
occurs  on  the  crystals  as  the  water  vapor  in  the  solution 
diffuses  to  the  crystal  and  freezes.  The  method  described  above 
ensures  that  no  water  vapor  enters  the  solution,  and  thus  clean 
formvar  slides  with  clear  ice  crystal  replicas  can  be  produced 
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every  time.  The  ice  crystal  and  droplet  concentrations  were 
calculated  using  the  equations  described  in  Appendix  I. 

Tne  charge  transfer  experiments  were  performed  by  moving 
targets  of  length  14  an  through  a  mixed  cloud  of  water  droplets 
and  ice  crystals.  In  the  large  ice  crystal  experiments  the 
target  diameters  were  5  mm  and  0.5  mm.  The  targets  were  mounted 
vertically  on  a  frame  attached  to  a  central  shaft  as  shown  in 
Figure  2.  The  shaft  could  be  rotated  so  that  the  targets  moved 
with  their  axes  normal  to  the  direction  of  motion  and  riming 
occurred  on  the  leading  edge  of  the  targets.  The  rod  speeds  used 
were  in  the  range  2-4  m  s  ^  but  some  experiments  were  performed 
at  a  speed  of  10  m  s  ^ .  Shielding  cups,  6  an  high,  were  placed 
at  both  ends  of  the  target  to  prevent  a  rime  bridge  from  forming 
and  linking  the  targets  to  earth.  The  shielding  cups,  Figure  2, 
were  grounded  to  remove  any  charges  acquired  during  an  experimen¬ 
tal  run. 

The  targets  were  connected  to  earth  via  a  sensitive  charge 
amplifier.  The  circuit  is  shown  in  Figure  3.  The  amplifier 
power  and  output  leads  were  taken  out  through  slip  rings 
positioned  at  the  bottom  of  the  shaft  under  the  cloud  chamber. 
The  555  timer  converts  the  positive  supply  voltage  to  a  negative 
supply  voltage  which  is  regulated  by  the  79105  voltage  regulator. 
The  charge  amplifier  is  a  CA3140S  chip  which  will  operate  down  to 
temperatures  of  -55°C.  Both  positive  and  negative  supply  rails 
are  required  so  that  the  amplifier  can  respond  to  an  input 
current  of  either  polarity.  The  background  noise  level  of  the 
charge  amplifiers  was  about  1  mV  so  the  charge  acquired  by  the 
target  was  measurable  if  the  charging  rate  was  greater  than  0.1 
pC  s  1.  The  supply  voltage  was  fused  at  5  volts.  The  amplifier 
has  a  sensitivity  of  1  mV  =  10  ^  A;  charge  flowed  to  earth  with 


8 


Figure 


a  time  constant  of  1  second.  Chart  recorders  were  used  to  record 
the  charging  currents  from  the  two  targets.  In  addition,  the 
rotating  frame  carried  two  other  probes  for  measuring  the  liquid 
water  content  of  the  cloud  and  the  rime  accretion  rate  of  a  5  mm 
target . 

The  rime  accretion  rate,  RAR,  is  the  amount  of  rime  collec¬ 
ted  by  the  target  per  unit  time  per  unit  area,  Jarayatne  and 
Saunders  (1985),  and  was  found  to  be  an  important  parameter  in 
controlling  the  charge  transfer.  As  the  droplets  freeze,  the 
rime  surface  is  warmed  due  to  the  release  of  latent  heat.  This 
rise  in  temperature  above  ambient  was  measured  with  an  accuracy 
of  0.13°C  by  a  platinum  resistance  temperature  sensor  of  diameter 
5  mm  which  thus  experienced  the  same  liquid  water  content  as  the 
5  mm  target.  By  using  the  equations  of  Ludlam  (1951)  and  Macklin 
and  Payne  (1967,  which  are  described  in  Appendix  II,  the  effec¬ 
tive  liquid  water  content,  Ew,  which  is  made  up  from  the  cloud 
droplets  large  enough  to  hit  the  target,  may  be  calculated  from 
the  temperature  elevation  of  the  rime.  Hence  the  rime  accretion 
rate  can  be  determined.  There  is  no  error  in  the  measurement 
introduced  by  the  presence  of  ice  crystals  in  the  cloud,  as  they 
do  not  release  latent  heat  on  striking  the  target.  The  rotating 
frame  also  carried  a  heated  wire  which  was  used  to  measure  the 
total  liquid  water  content  of  the  cloud. 

In  conjunction  with  the  rotating  targets,  there  was  also  a 
stationary  target  connected  to  a  charge  amplifier.  The  target, 
of  length  2  an,  was  situated  inside  a  tube  of  diameter  2  an  which 
protruded  a  short  distance  into  the  cloud  chamber.  A  regenera¬ 
tive  pump  drew  the  cloud  past  the  target  at  speeds  of  up  to  110  m 
s  which  is  typical  of  some  aircraft  speeds.  Three  separate 
sets  of  experiments  were  performed  with  a  clean  metal  target,  a 


rimed  target  and  a  riming  target. 

Thus  during  each  experimental  run  the  charging  current  to 
two  targets,  the  rime  accretion  rate  and  the  liquid  water  content 
of  the  cloud  could  be  measured  simultaneously. 
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SECTION  3 


A  NEW  HOT-WIRE  FOR  MEASURING  CLOUD  LIQUID  WATER  CONTENT 

Introduction 

The  measurement  of  wind  velocity  by  monitoring  the  power 
dissipation  in  a  hot-wire  was  first  suggested  by  Kennelly  (1909) 
and  tested  experimentally  by  Morris  (1912).  King  (1914)  detailed 
extensively  the  convection  of  heat  from  small  cylinders  in  a 
fluid  flow,  using  Boussinesq's  (1905)  equation  for  heat  loss  to  a 
stream  of  fluid.  He  suggested  that  there  were  applications  for 
'hot-wire  anemometry'. 

Merceret  and  Schricker  (1975),  from  the  work  of  Wynngaard 
and  Lumley  (1967),  adapted  a  hot-wire  anemometer  for  measuring 
the  liquid  water  content,  l.w.c,  of  a  cloud.  They  found  that  the 
constant  temperature  hot-wire  nimbiometer  was  far  superior  to  the 
constant  current  hot-wire  devices  such  as  the  Johnson-Will iams 
instruments.  The  J-W  instruments  will  only  respond  accurately  to 
droplets  <30  pm  and  have  a  high  power  consumption,  difficult 
calibration  and  a  time  constant  longer  than  one  second.  Merceret 
and  Schricker  found  a  relationship  between  the  power  used  and  the 
mass  of  water  intercepted  by  the  wire  with  time.  With  their 
analysis  they  did  not  need  to  supply  heat  of  vaporisation  as  the 
droplets  were  not  evaporated  due  to  their  low  operating  tempera¬ 
ture.  They  calculated  the  proportion  of  wire  that  is  wet  at  any 
instant  in  time,  and  so  only  needed  to  ’cut’  the  droplets  rather 
than  evap-orate  them,  and  hence  there  will  be  little  spectral 
filtering  over  a  wide  range  of  droplet  sizes. 

King,  Parkin  and  Handsworth  (1978)  also  preferred  the 
constant  temperature  device  for  use  on  an  aircraft.  They  used  a 
coil  of  1.7  mm  diameter  and  length  34  mm.  They  found  that  slave 
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coils  on  either  end  were  necessary  to  maintain  a  fairly  constant 
temperature  along  the  master  coil.  The  slave  coils  are  17  mm  in 
length  and  1.7  mm  in  diameter  and  do  not  form  part  of  the  sensing 
probe.  They  used  a  simple  control  circuit  to  maintain  a  constant 
temperature.  An  analysis  of  a  similar  control  circuit  was  pre¬ 
sented  by  Freymuth  (1977).  Various  operating  temperatures  for 
the  wire  were  tried,  ranging  from  76°C  to  178°C  in  both  wind 
tunnel  and  cold-room.  An  operating  temperature  in  the  region  of 
80-90°C  was  recommended  for  two  reasons: 

(a)  At  some  altitudes  the  boiling  point  of  water  can  be  as  low 
as  90 °C. 

(b)  If  the  temperature  is  in  excess  of  100°C  then  evaporation 
can  take  longer  due  to  the  formation  of  an  insulating  layer 
between  the  droplets  and  the  wire. 

Later,  however.  King  et  al .  (1981)  state  that  the  wire 

temperature  can  be  as  high  as  160°C  before  any  nucleate  boiling 
of  the  drops  is  observed.  This  is  due  to  a  30  m  thick  covering 
of  epoxy  resin  on  the  wire  which  has  a  surface  temperature  of 
75°C  for  a  wire  temperature  of  160°C.  King  et  al .  (1984)  state 

that  with  an  epoxy  covered  wire  the  wire  had  to  be  at  least  60°C 
warmer  than  the  unshielded  wire  in  order  to  achieve  the  same 
response.  They  also  state  that  the  optimum  operating  temperature 
of  the  unshielded  probe  is  about  180-200°C,  an i  therefore  the 
shielded  version  would  operate  at  about  2408C,  but  this  leads  to 
failure  of  the  device.  At  a  temperature  of  2208C  the  response  of 
the  shielded  probe  is  about  0.9  that  of  the  unshielded  probe. 
Bradley  and  King  (1980  )  and  King  et  al  .  (1981,  1984  )  showed  the 

frequency  response  of  the  electronic  circuit  and  presented  de¬ 
tails  of  extensive  wind  tunnel  and  cold-room  tests  performed  on 


12 


the  probe . 

The  following  sections  outline  the  principles  behind  the 
operation  of  a  hot-wire  device  and  the  calibration  of  such  a 
device . 

General  Operating  Principles  of  a  Hot-Wire  Device 

If  a  hot  object  is  placed  in  a  colder  fluid  flow,  then  it 
will  lose  heat  until  it  is  in  temperature  equilibrium  with  the 
fluid.  If,  however,  a  current  is  supplied  to  the  object,  in  this 
case  a  metal  wire,  then  a  constant  temperature  can  be  maintained 
and  the  power  supplied  gives  an  indication  of  the  fluid  velocity 
in  the  case  of  hot-wire  anemometry,  or  cloud  liquid  water  con¬ 
tent,  l.w.c.,  in  the  case  of  hot-wire  nimbiometry. 

King's  law,  after  L  V  King  (1914),  states  that  for  small 
cylinders  the  general  heat  transfer  law  is: 

P  *  (A  +  BUN>  <Tw-Ta) 

where:  P  is  the  rate  of  cooling  of  the  object  and  has  the 

dimensions  of  power,  A, B , N  are  the  'King's  law  constants',  U  is 
the  velocity  of  the  fluid,  T  ,T  are  the  wire  and  air  tempera- 
tures  respectively. 

King  et  al.  (1978),  however,  state  that  the  power  necessary 
to  maintain  the  probe  wire  at  a  particular  temperature  can  be 
determined  by  the  dimensions  of  the  probe,  and  is  given  by: 

P  =  ldvw[L+c(Tw-Ta )]  +  lk(Tw-Ta)Nu 
where:  1  is  the  wire  length,  d  is  the  wire  diameter,  v  is  the  air 
velocity,  w  is  the  liquid  water  content,  L,  c  are  the  latent  and 
specific  heat  of  water  respectively,  Tw'Ta  are  the  wire  and  air 
temperatures  respectively,  k  is  the  thermal  conductivity  of  air, 
and  Nu  is  the  Nusselt  number  for  the  heat  transferred  from  the 
wire.  The  first  term  is  known  as  the  'wet'  term,  and  is  the  heat 
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required  to  warm  the  supercooled  water  droplets  from  temperature 
to  and  then  to  evaporate  them.  The  second  or  'dry'  term 
represents  the  heat  transferred  to  the  cooler  air  flowing  over 
the  wire. 

There  were  two  possible  approaches  to  the  proble.  of 
positioning  the  probes:  (a)  a  rotating  probe  or  (b)  a  stationary 
probe  with  the  cloud  being  drawn  past  it.  The  nature  and  cali¬ 
bration  of  both  probes  and  the  electronic  control  circuitry  are 
detailed  in  the  following  sections. 

The  Electronic  Control  Circuit 

The  novel  control  circuit,  shown  in  Figure  4,  uses  a 
switching  signal  to  regulate  the  operating  temperature,  whereas 
other  workers  have  used  a  linear  circuit  to  perform  this 
function.  The  wire  is  heated  by  applying  to  it  a  chain  of  fixed 
voltage  pulses  whose  mark-space  ratio  is  varied  to  control  the 
mean  electrical  power  dissipation  in  the  wire.  As  the  loading  on 
the  wire  increases  due  to  an  increase  in  cloud  liquid  water 
content,  or  a  decrease  in  air  temperature,  the  mark-space  ratio 
increases . 

The  555  time,  IC1,  is  running  as  an  astable  multivibrator 
which  produces  a  square  wave  output  with  a  high  time  of  0 . 7  ms . 
The  4027,  IC2,  is  a  J-K  Flip-Flop  which  changes  from  a  low  state 
output  of  0  V  to  a  high  state  output  of  12  V  on  pin  2.  When  pin 
4  is  high,  pin  2  is  high,  and  vice-versa.  The  4027  triggers  the 
transistor  which  in  turn  applies  the  rail  voltage  of  12  V  across 
the  bridge  formed  by  R?,  Rg,  R1Q  and  the  wire  (Rg).  The  wire  is 
heated  until  the  bridge  reaches  a  balance.  The  voltage  compara¬ 
tor,  IC3,  monitors  the  voltage  across  the  two  arms  of  the  bridge 
and,  when  balance  is  reached,  it  resets  IC2  which  turns  TRl  off. 
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The  electronic  control  circuit 


This  sequence  occurs  for  every  cycle  of  the  555  timer.  The 
frequency  is  set  sufficiently  high  so  that  the  cycle  time  is 
short  compared  with  the  thermal  response  time  of  the  wire.  This 
short  cycle  time  ensures  that  the  fall  in  the  wire  temperature  is 
small,  compared  to  its  operating  temperature  during  the  off  time 
of  the  transistor. 

The  coupling  capacitors  across  pins  1  and  8  (555),  pins  1 
and  8,  and  1  and  4  (311)  are  simply  to  smooth  out  any 
irregularities  in  the  power  supply.  When  the  circuit  is  initial¬ 
ly  switched  on  there  needs  to  be  an  offset  of  a  few  millivolts  on 
pin  3  of  the  311  to  send  the  311  output  high  and  subsequently 
turn  the  4027  on. 

The  transistor  TR1 ,  TIP  146,  is  a  pnp  type  transistor.  A 
pnp  is  used  in  preference  to  a  npn  type  due  to  the  inverted 
signal  from  pin  2  of  the  4027,  and  so  when  the  4027  switches  off 
there  is  less  heat  dissipation  in  the  transistor  due  to  the 
reversal  in  base  current  when  the  4027  switches  off. 

The  voltage  across  R8  (wire)  is  a  square  wave  signal;  the 
'on'  time  increasing  as  air  or  water  droplets  are  blown  across 
the  wire.  The  actual  output  is  at  point  A,  after  the  filter 
circuit  which  effectively  removes  the  square  wave.  There  is  a 
voltage  noise  on  the  signal  of  approximately  15  mV,  which  is 
primarily  due  to  noise  from  the  slip-rings. 

To  set  the  operating  temperature  of  the  wire,  two  poten¬ 
tiometers  are  used;  the  first  to  turn  on  the  4027  and  heat  the 
wire  to  almost  the  correct  temperature,  and  the  second,  Ik,  is 
used  to  set  accurately  the  temperature  due  to  its  greater  sensi¬ 
tivity.  At  the  required  operating  temperature,  the  wire  resis¬ 
tance  is  Rw.  The  voltage  applied  across  the  wire,  V  ,  during  the 
on  time  of  the  transistor  is; 
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V  =  R  V  /{R  +R_ )  (  1 ) 

w  w  s  w  7 

where  Vg  is  the  supply  voltage 

The  power  dissipated  in  the  wire  during  this  time  is: 

Pw  =  Vw2/Rw  =  RwVs2/(Rw+R7)2  (2) 

If  the  on  time  is  of  duration  t  and  the  total  cycle  is  t,  then 

on 

the  average  power,  p  ,  dissipated  in  the  wire  is: 

Fw  “  Von^  *  RwVs2/(Rw  +  R7)2  <3> 

The  mean  voltage  across  the  wire,  v^,  is: 


/  =  V_  t  /t  =  R^  Ve  t  / <R  +R, )t  (4) 

w  w  on  w  s  on  w  7 


In  practice  V  ^is  derived  via  a  low  pass  filter  at  point  A.  Thus 
from  equations  3  and  4 

VVsVVR7> 

Thus  there  is  a  linear  relationship  between  and  the  power 
dissipation.  Therefore,  the  fractional  changes  in  Vw  are  equal 
to  the  fractional  changes  in  power  dissipation,  that  is; 

dVVw  =  dVPw  (6) 

In  the  conventional  circuits  used  by  King  and  others,  the  power 


dissipation  is  given  by: 


therefore : 


P  =  V  VR 
w  w  w 


dVw/Vw  =  °-5  dpw/pw  <8) 

From  equations  6  and  8  it  can  be  seen  that  the  pulse  width 

modulated  control  provides  an  output  which  has  twice  the  sensiti¬ 
vity  to  power  dissipation  changes  of  a  conventional  circuit. 
Also,  there  is  a  reduction  of  dissipation  in  the  controlling 
transistor  which  leads  to  reduced  overall  power  and  instrument 
weight.  As  indicated  by  King,  conventional  circuits  require 
adjustment  of  the  frequency  compensation  to  optimise  performance 
and  avoid  instability.  This  is  not  required  with  pulse  width 
modulated  control. 
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In  order  to  maintain  a  constant  temperature  across  the 
probe,  King  et  al.(1978)  had  to  use  slave  coils  on  both  ends  of 
the  probe.  The  coils  were  not  part  of  the  sensing  circuit. 
However,  in  this  case,  the  wire  is  of  a  very  small  diameter  and 
of  relatively  short  length,  and  so  the  temperature  gradient 
losses  that  King  encountered  can  be  ignored. 

Calibration  and  Results 

The  probe  was  operated  at  a  wire  temperature  of  76°C.  The 
reason  being  that  at  higher  wire  temperatures  there  occurs  film 
boiling  at  the  surface  of  the  wire  causing  an  insulating  layer  of 
vapor  between  the  droplet  and  the  wire.  The  responses  at  76°C 
and  98  °  C  are  7.27  g  m  ^  V  ^  and  9.33  g  m  ^  V  ^  respectively, 
indicating  that  there  is  a  better  response  at  the  lower  operating 
temperature.  At  still  lower  temperatures,  56 °C  and  36°C,  there 
was  evidence  of  the  wire  becoming  saturated  with  water  droplets 
leading  to  a  longer  evaporation  time.  Located  on  the  opposite 
arm  to  the  probe  was  a  collecting  rod  of  length  100  mm  and 
diameter  1.5875  mm.  The  LWC  of  a  cloud  was  deduced  from  the 
weight  of  rime  accreted  on  the  rod  after  it  had  been  rotated  for 
a  time  ( Mackl in , 1962  ) .  The  LWC  values  achieved  using  this  method 
were  corrected  for  the  collision  efficiency  of  the  rod.  The 
droplet  spectra  (obtained  with  a  Knollenberg  optical  probe)for 
two  different  clouds  were  used  to  average  the  collision  efficien¬ 
cy  (CE)  of  the  rod  over  all  droplet  sizes  by  using  the  equation: 

CE  =  (SNdEd,/Ntot  (9) 

where  Ntot  is  the  total  number  of  droplets/am^ ,  is  the  number 

of  droplets  in  size  range  Ad  and  is  the  collision  efficiency 
of  the  rod  for  droplets  of  diameter  d  yin  (Ranz  and  Wong,  1952). 

By  using  this  analysis  the  mean  efficiency  of  the  rod,  57%, 


17 


agreed  within  1%  for  the  two  different  clouds.  All  of  the  clouds 
were  produced  by  a  nebuliser,  and  they  had  a  broad  spectrum  of 
droplet  sizes  with  a  mean  droplet  diameter  of  12  pm. 

The  initial  tests  were  performed  at  a  temperature  of  -9°C 
and  then  repeated  at  -3°C  and  14°C.  These  results  are  shown  in 
Figure  5,  and  it  can  be  seen  that  the  instrument  gives  a  linear 
response  dependent  on  the  LWC  of  the  cloud  and  its  ambient  air 
temperature.  The  lines  have  an  equation: 

LWC  =  7.27V  +  0.1575T.  (10) 

As  the  lines  are  parallel,  it  is  possible  to  use  just  one  line  and 
introduce  a  temperature  factor  such  that 

V  .  =  AV  +  0 . 02 2 AT  (  11  ) 

C3  1 

where  V  ,  is  the  voltage  on  the  calibration  line,  AV  is  the 
recorded  output  above  1.5V  and  ^T  is  the  temperature  difference 
from  -9°C.  (Warmer  than  -9°C  is  positive  and  colder  than  -9°C  is 
negative ) . 

The  calibration  remained  unchanged  if  a  new  wire  of  the  same 
specification  and  length  was  used.  The  calibration  could  be 
checked  very  easily  by  rotating  the  probe  in  clean  air  at  a  known 
temperature  while  monitoring  the  output  and  comparing  this  with 
the  original  calibration. 

Figure  6  shows  the  clear  air  output,  V  against  the  output 
Vc,  obtained  when  the  probe  is  rotated  in  a  crystal  cloud.  It 
was  important  that  the  cloud  consisted  purely  of  ice  crystals,  as 
the  presence  of  a  few  droplets,  even  10  at  can  give  rise  to 
considerable  errors.  This  experiment  was  repeated  at  different 
air  temperatures  and  with  different  sized  crystals  which  were 
obtained  by  different  methods  of  seeding  a  supercooled  cloud  by 
introducing  into  it  a  thin  rod  which  had  been  dipped  in  liquid 
nitrogen.  Small  crystals  were  obtained  by  leaving  the  rod  in  the 
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cloud  for  a  relatively  long  time  of  around  4  s.  In  the  majority 
of  these  experiments  the  cloud  output,  Vc,  equalled  the  clear  air 
output,  V  ,  giving  a  1:1  correlation.  In  the  worst  case,  there 

a 

was  an  output  of  5  mV  above  the  value  of  V  which  corresponds  to 
a  LWC  of  0.01  g  m  ^  which  is  well  within  the  accuracy  of  the 
device.  This  means  that  in  a  mixed  cloud  of  ice  crystals  and 
supercooled  water  droplets  the  crystals  will  be  ignored  and  the 
recorded  output  will  be  a  true  measure  of  the  LWC,  with  no  error 
introduced  by  the  presence  of  the  crystals  irrespective  of  the 
crystal  diameter. 

Concl usion 

This  hot-wire  device  is  suitable  for  use  in  the  laboratory, 
as  it  needs  only  a  simple  calibration  and  the  data  may  be  presen¬ 
ted  in  the  form  of  a  voltage-time  graph.  The  probe  itself  is 
robust,  easy  to  construct  and  inexpensive.  The  total  error  was 
determined  by  taking  the  square  root  of  the  sum  of  the  square  of 
the  individual  errors.  The  source  of  these  errors  were:  the 
mass  of  rime  collected,  the  diameter  and  velocity  of  the  collec¬ 
ting  rod  and  the  collection  time.  The  most  significant  error  was 
in  the  mass  of  accreted  rime.  The  instrument  has  a  response  time 
of  less  than  0.5  s  and  is  accurate  to  4%  at  l.Og  m  ^ ,  4.6%  at 
0.5  g  m  ^  and  13.2%  at  0.06  g  m~^. 
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SECTION  4 


THE  COLLECTION  EFFICIENCY  OF  A  RIME  COVERED  TARGET 

Introduction 

It  is  the  purpose  of  this  part  of  the  study  to  investigate 
whether  soft-hailstones  have  an  increased  collection  efficiency 
for  water  droplets  due  to  the  presence  of  rime  feathers.  The 
study  also  investigates  the  dependence  of  the  collection  effi¬ 
ciency  on  temperature  and  impact  velocity. 

The  accretion  of  supercooled  water  droplets  upon  cloud  par¬ 
ticles  is  known  as  riming,  and  is  the  dominant  growth  process  for 
large  hydrometeors.  It  is  important  to  know  the  values  of  the 
collection  efficiency  of  growing  soft-hailstones  for  supercooled 
water  droplets  through  which  they  fall,  because  the  charge 
transfer  due  to  ice  crystal  collisions  appears  to  be  affected  by 
the  presence  of  liquid  water  on  the  riming  target. 

Ranz  and  Wong  (1952)  modelled  theoretically  the  collision 
efficiency  of  a  smooch  cylindrical  target  for  smoke  particles  of 
a  known  diameter.  Many  workers,  including  Jayaratne  et  al . 
(1983),  have  used  their  results  to  determine  the  growth  rate  of  a 
soft-hailstone  when  it  moves  through  a  cloud  of  supercooled 
water  droplets  of  known  drop  size  distribution  and  number  concen¬ 
tration.  The  collision  efficiency  is  assumed  to  be  the  same  as 
the  collection  efficiency  for  this  case,  because  the  droplets 
freeze  and  stick  on  impact  with  the  collector.  Jayaratne  et  al . 
simulated  thunderstorm  conditions  in  the  laboratory,  and  drew  ice 
crystals  and  supercooled  droplets  past  a  rimed  target  which 
represented  a  soft-hailstone  of  5mm  diameter  falling  at  3  m  s 
for  this  case,  according  to  Ranz  and  Wong,  the  target  does  not 
collect  any  droplets  small  than  6  yum  diameter. 


The  liquid  water  content,  w,  of  a  supercooled  cloud  can  be 
measured  from  the  amount  of  rime  which  is  accreted  on  a  moving 
target  in  a  given  time,  as  shown  by  Brun  et  al.  (1955),  Rogers  et 
al .  (1983  ),  Saunders  et  al.(1984  ),  and  in  this  report.  The 

method  is  only  possible  if  the  drop  size  distribution  of  the 
cloud  and,  also,  the  collection  efficiency  of  the  target  for  the 
droplets  is  known  accurately.  The  density  of  rime  grown  in  this 
way  was  investigated  by  Macklin  (1962),  and  Saunders  and  Zhang 
(1987).  In  both  cases,  the  targets  were  exposed  to  a  stream  of 
water  droplets  which  froze  on  impact  and  formed  the  rime.  From 
the  geometrical  dimensions  and  weight,  the  rime  density  was 
calculated.  This  method  provided  a  measure  of  the  effective 
liquid  water  content  in  the  cloud,  Ew,  where  E  is  the  collection 
efficiency.  Measurements  of  the  density  of  rime  are  only  con¬ 
cerned  with  the  amount  of  water  that  collides  with  the  target, 
and  so  this  method  removes  the  need  to  calculate  E  accurately. 
The  results  of  Saunders  and  Zhang  compared  favorably  with  those 
of  Macklin.  At  low  impact  speeds  and  temperatures,  the  growth  of 
rime  'feathers'  is  enhanced.  The  term  'feathers'  was  used  by 
Macklin  to  describe  the  rime  appearance  following  the  freezing  of 
individual  droplets  on  top  of  one  another,  so  that  the  rime 
deposit  grew  out  from  the  ice  surface.  Under  these  conditions,  a 
riming  target  will  not  be  of  uniform  cross  section,  and  the 
collision  efficiency  theory  of  Ranz  and  Wong  may  not  be  strictly 
appl ica  ble . 

Experimental  Method 

The  experiments  were  performed  in  tr  ud  chamber  situated 

within  the  large  cold-room.  The  aim  c  the  experiments  was  to 
determine  whether  small  droplets  are  col.  cted  by  rime  feathers. 
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In  order  to  achieve  this,  the  droplets  were  formed  from  a  strong 
salt  solution  and  then  drawn  past  a  previously  rimed  target  of 
length  53mm.  The  target  diameter  was  5mm,  which  is  the  same  size 
as  that  used  in  the  charging  experments.  A  plastic  target  was 
used  in  preference  to  a  metal  one,  because  rime  scraped  from  a 
metal  target  contained  tiny  metal  shavings  which  subsequently 
contaminated  the  rime  solution.  The  target  material  did  not 
affect  the  nature  of  the  riming. 

The  presence  of  the  salt  droplets  in  the  rime  was  detected 
by  a  Kent  conductivity  meter.  If  the  droplets  hit  the  rime 
feathers,  then  the  pre-existing  pure  rime  became  contaminated 
with  the  salt  solution.  As  contamination  increased  due  to 
increasing  exposure  to  the  salt  droplet  cloud  or  increasing 
collection  efficiency  with  feather  development,  the  rime 
conductivity  also  increased.  The  conductivity  meter  was  first 
calibrated  with  known  concentrations  of  NaCl  solution.  The 
results,  presented  in  Figure  7,  show  that  the  meter  has  a  linear 
response  with  increasing  salt  solution  concentration.  This 
linearity  continues  up  to  100  g  1  ^ ,  at  which  point  the  solution 
starts  to  become  supersaturated  with  salt.  The  resonse  of  the 
instrument  was  0.05  pS  cm^  ,  which  corresponds  to  a  change  of 
0.006  pi  in  the  amount  collected.  The  salt  solutions  were  made 
up  with  NaCl  and  deionised  distilled  water,  so  as  to  minimise  the 
effect  of  any  impurities,  and  to  ensure  that  the  solutions  were 
of  a  consistent  conductivity.  New  standard  salt  solutions  were 
used  every  day  because  the  concentration  of  the  solution 
increased  due  to  evaporation  of  the  water.  The  concentration 
used  was  10  gl  ”1  a  concentration  of  1  g  1  _1was  tried  but  the 
contamination  of  the  rime  was  below  the  threshold  of  detectabili¬ 
ty  of  the  conductivity  meter.  The  meter  was  checked  periodical- 
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introduced  water  vapor  into  the  cloud  chamber  at  a  constant  rate 
to  produce  droplets  which  subsequently  cooled  to  the  ambient 
temperature.  The  cloud  had  a  modal  drop  size  of  12  pm  and  the 
liquid  water  content  was  typically  between  1-2  g  m  .  The  ex¬ 
periments  were  performed  at  temperatures  in  the  range  -5  to  -20°C 
with  impact  speeds  of  3  and  5  m  s-^.  The  cloud  was  drawn  past 
the  target  at  a  constant  rate  for  5  minutes,  which  was  found  to 
be  sufficient  time  for  the  rime  feathers  to  grow.  The  rimed 
target  was  then  exposed  to  the  salt  droplet  cloud  produced  by  a 
snail  Unicorn  nebuliser.  This  cloud  was  drawn  past  the  target 
for  5  minutes,  at  the  same  velocity  as  that  used  in  the  pure 
riming  part  of  the  experiment.  All  the  contaminated  rime  was 
then  removed  from  the  target  and  diluted  with  20  ml  of  deionised 
distilled  water . 

The  drop  size  distribution  of  the  salt  cloud  was  analysed  by 
a  Knollenberg  Forward  Scattering  Spectrometer  Probe  (FSSP).  The 
modal  diameter  was  about  10  pm,  and  a  typical  drop  size  distribu¬ 
tion  is  shown  in  Figure  8.  As  can  be  seen,  there  were  some 
droplets  larger  than  20  pm  which  were  not  desirable  but  were 
unavoidable.  The  presence  of  these  droplets  was  accounted  for  in 
the  analysis  of  the  data. 

The  theoretical  collection  efficiencies,  shown  in  Figure  9, 
were  calculated  from  the  equations  of  Loffler  and  Muhr  (1972), 
for  a  smooth  collector  of  diameter  5mm,  cloud  temperature  of 
-10®C,  with  air  velocities  of  3  and  5ms1.  The  Ranz  and  Wong 
efficiencies  used  by  Jayaratne  et  al.  over-estimated  the  collec¬ 
tion  efficiency  for  small  droplets  by  about  10%,  compared  with 
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the  Loffler  and  Muhr  equations.  The  collection  efficiency  in¬ 
creases  slightly  with  decreasing  temperature,  due  to  the  lower 
viscosity  of  air  at  lower  temperatures.  For  the  3  m  s"1  air 
velocity,  it  can  be  seen  that  the  collection  efficiency  for 
droplets  of  diameter  between  5  and  8  pm  is  less  than  1.0%.  These 
theoretical  collection  efficiencies  are  for  smooth  collectors, 
and  do  not  take  into  account  any  turbulence  due  to  the  roughness 
of  the  surface. 

The  total  volume  of  salt  water  droplets  collected  by  the 
sof t -hai lstone  target  can  be  calculated  from  the  following 
formula  . 

V.  =  (Xr)(V  +  V  ,  +  20xl0-3 ) / ( X . , )  (1) 

a  Led  N 

where  X..  is  the  concentration  of  salt  solution  in  the  nebuliser, 
h 

(g/1),  XL  is  the  concentration  of  salt  solution  in  the  rime  plus 
20ml  of  water,  (g/1);  is  the  total  volume  of  salt  solution 

droplets  which  collide  with  the  target,  (1);  Vc  is  the  volume  of 
pure  rime  on  the  target,  (1), 

Now,  (V  +  V.)  is  much  less  than  20ml,  so  from  (1) 
c  a 

Vd  =  20x10 ”3  X  /(X^  (  2) 

Thus  from  the  salt  concentration  of  the  rime  liquid  and  the 
initial  salt  solution,  the  volume  of  salt  cloud  which  hits  the 
rime  target  can  be  calculated.  From  the  Lcffler  and  Muhr 
formula  and  the  measured  droplet  sire  distribution,  the  total 
volume  of  salt  droplets  that  can  collide  with  the  target  is  VD> 
So,  if  >  Vp  then  the  difference  is  made  up  from  the  extra 

droplets  that  collide  which  are  not  accounted  for  by  the  collec¬ 
tion  efficiency  formula  for  smooth  targets. 
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Figure  3 
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Results 


The  experiments  were  performed  at  different  temperatures 
between  -5  and  -20°C  at  a  constant  air  velocity  of  3  m  s-^.  The 
pure  droplet  cloud  was  drawn  past  the  target  for  5  minutes  before 
being  exposed  to  the  salt  droplets  for  5  minutes.  The  validity 
of  the  equations  described  in  the  previous  section  was  confirmed 
by  varying  the  concentra t ion ,  X^,  of  the  nebuliser  solution. 
Figure  10  shows  that  the  volume  collected  is  approximately  the 
same  for  all  values  of  X^.  The  only  limits  on  the  solution 
concentration  were  of  detectability  at  low  values  and  of  supersa¬ 
turation  at  high  values. 

Figure  11  shows  the  volume  of  salt  droplets  collected  by  the 
rime  against  cloud  temperature.  At  -5°C  the  volume  collected  is 
around  28.9  pi  increasing  to  37.1  pi  at  a  temperature  of  -20°C. 

The  time  of  exposure  of  the  target  to  the  pure  droplet  cloud 
was  varied  from  2  to  7  minutes  at  a  temperature  of  -10° C.  The 
subsequent  increase  in  the  collection  efficiency  of  the  target 
for  the  salt  droplets  is  reflected  in  the  volume  of  the  salt 
droplets  collected  in  5  minutes.  Figure  12.  As  the  exposure  time 
increased,  the  number  and  size  of  the  rime  feathers  also  in¬ 
creased  for  the  first  6  minutes  of  collection.  Beyond  this  time 
the  collection  efficiency  remains  approximately  constant,  indica¬ 
ting  no  further  increase  in  feather  concentration. 

In  further  experiments,  performed  under  '.he  same  conditions 
as  above,  the  pure  droplet  cloud  was  drawn  past  the  target  at  a 
velocity  of  5  m  s  ^  for  5  minutes.  The  rime  acquired  the  kernel 
appearance  described  by  Macklin  (1962)  but  the  rime  feathers  were 
not  as  visible  as  in  the  above  case.  Figure  13  shows  the  salt 
droplet  volume  collected  in  5  minutes  against  temperature.  At 
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this  higher  impact  velocity,  the  increase  in  the  amount  collected 
with  decreasing  temperature  is  less  than  at  3  m  s-^.  The  reasons 
for  this  will  be  discussed  in  the  following  section. 

Discussion 

It  has  been  hypothesied  by  previous  workers,  Macklin  (1962), 
and  Dong  and  Hallett  (1986),  that  as  a  hailstone  grows  due  to  the 
accretion  of  water  droplets,  the  collection  efficiency  for  those 
droplets  increases  due  to  the  formation  of  rime  feathers. 

Mossop  (1984)  showed  that  in  the  production  of  laboratory 
clouds  the  drop  size  distribution  is  affected  by  the  size  of  the 
water  vapor  inlet  nozzle.  It  can  be  seen  from  Figure  14a  that  as 
the  nozzle  size  decreases  the  drop  size  distribution  shifts  to  a 
higher  concentra tion  of  smaller  droplets.  This  was  confirmed  by 
Jayaratne  and  Saunders  (1985),  Figure  14b,  when  they  showed  that 
the  charging  of  a  hailstone  due  to  the  collision  and  separation 
of  ice  crystals  is  dramatically  affected  by  the  drop  size  distri¬ 
bution  present  in  the  cloud.  Figure  15a  shows  the  charging  which 
occurs  with  a  large  nozzle  of  diameter  32  mm  at  -10°,  whereas 
Figure  15b  shows  the  charging  under  similar  conditions  but  with 
the  small  droplet  cloud  produced  by  the  small  nozzle  of  diameter 
5  mm.  In  the  second  case,  when  negative  charging  occurs,  the  ice 
crystals  have  grown  at  the  expense  of  the  droplets  which  have  a 
maximum  size  of  6pm,  and  so  are  too  small  to  strike  the  target 
according  to  the  theory  of  Ran z  and  Wong.  Jayaratne  and  Saunders 
(1985),  and  Baker  et  al.(1987),  concluded  from  these  experiments 
that  it  is  the  surface  properties  of  the  interactants  which  are 
important  in  controlling  the  sign  of  the  charge  transfer. 
Charging  is  positive  if  the  hailstone  is  growing  faster  than  the 
ice  crystals,  and  negative  if  growing  more  slowly.  The  results 
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of  the  present  study  show  that  some  riming  due  to  the  small 
droplets  may  occur  due  to  the  presence  of  rime  feathers  on  the 
target,  and  this  brings  the  charge  sign  hypothesis  into  question. 
It  can  be  shown,  however,  by  using  the  formulae  of  Baker  et  al., 
that  the  amount  of  water  accreted  by  the  hailstone  is  insuffi¬ 
cient  to  increase  the  hailstone  growth  rate  above  that  of  the  ice 
crystals,  and  so  the  sign  of  the  charging  is  unaffected.  A 
detailed  account  of  the  relative  growth  rates  of  the  cloud  parti¬ 
cles  is  given  in  Appendix  III. 

The  appearance  of  the  rimed  surface  under  varying  conditions 
has  been  described  in  detail  by  Macklin  (1962),  Stallabrass 
(1978)  and  Lowaski  (1983).  At  low  temperatures  and  low  impact 
velocities,  the  droplets  freeze  without  significant  flattening 
which  enhances  the  growth  of  the  rime  feathers.  At  higher  tem¬ 
peratures,  the  rime  appears  as  smooth  clear  ice.  As  the  drplets 
freeze  they  release  latent  heat  of  fusion  which  warms  the  rime 
surface  above  the  ambient  cloud  temperature.  Using  the  results 
of  Ludlam  (1951),  it  can  be  shown  that  if  a  cloud  of  temperature 
-5°C  and  liquid  water  content  1-  g  m"*  is  drawn  past  a  5mm  diame¬ 
ter  target  at  5  m  s  *  then  the  surface  temperature  of  the  rime 
increases  to  -2.9°C.  Under  similar  conditions,  but  with  a  cloud 
temperature  of  -8°C,  the  rime  surface  temperature  will  be  -5.9“C. 
Thus,  conditions  of  wet  growth  are  not  achieved  in  these  experi¬ 
ments.  However,  the  droplet  freezing  time  is  temperature  depen¬ 
dent,  with  longer  freezing  times  at  higher  temperatures  giving 
time  for  the  freezing  droplets  to  flow  over  the  surface,  leading 
to  a  smooth  ice  surface.  Loffler  and  Muhr  assumed  that  the 
collector  is  smooth,  and  thus  their  theory  and  the  present  ex¬ 
perimental  results  agree,  with  a  near  zero  target  collection 
efficiency  at  temperatures  above  -8°C.  In  the  case  of  a  smooth 
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collector,  the  theoretical  collection  efficiency  increases  with 
increasing  air  velocity.  However,  if  the  collector  is  rough  then 
the  experiments  show  that  the  increase  in  the  volume  collected 
due  to  the  roughness  decreases  with  increasing  air  velocity,  as 
indicated  in  Figure  16.  The  volume  collected  at  -5°C  (3  m  ) 
and  -8°C  (5  m  s^  )  is  taken  as  zero,  and  the  relative  increase  in 
volume  collected  due  to  the  presence  of  the  rime  feathers  is 
shown  as  a  function  of  temperature.  It  can  be  seen  that  at  low 
speeds  where  feather  growth  is  enhanced,  the  volume  collected 
increases  substantial ly  at  lower  temperatures.  At  a  flow  speed 
of  5  m  s  *  the  feather  growth  is  inhibited,  and  thus  the  collec¬ 
tor  surface  is  essentially  smooth.  Macklin  noted  that  with  an 
increase  in  impact  velocity  the  rime  surface  appearance  changed 
from  a  'feathery'  to  a  'kernal'  type  due  to  gaps  in  the  rime 
being  filled-in  by  subsequent  droplets  impacting  at  the  higher 
velocity.  This  change  will  occur  at  lower  temperatures  for 
higher  values  of  liquid  water  content.  For  instance,  at  a  tem¬ 
perature  of  -150C  where  rime  feather  growth  is  prominent,  the 
surface  temperature  would  rise  to  -10°C  if  the  liquid  water 
content  was  increased  to  3  g  m  ^ .  The  growth  of  rime  feathers 
would  then  be  inhibited  and  the  collection  efficiency  of  the 
target  would  decrease. 

The  equivalent  collector  diameter  is  defined  as  that  smooth 
collector  diameter  whch  would  theoretically  (T,offler  and  Muhr) 
collect  the  same  volume  of  droplets  as  in  the  experiment.  At  low 
air  speeds  (3ms  ^),  where  feather  growth  occurs,  this  diameter 
is  3.2mm;  whereas,  at  5  m  ^  the  equivalent  collector  diameter 
is  4.8mm.  The  presence  of  rime  on  the  collector  increases  the 
collection  efficiency  for  all  droplet  sizes.  At  3m  s~*  the 
collection  efficiency  for  5  pm  droplets  has  increased  from  0.06% 
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for  a  smooth  5mm  collector  to  0.211  for  a  rough  surface  with  an 
equivalent  collector  diameter  of  3.2mm.  Table  1  shows  the  coll¬ 
ection  efficiency  for  a  5  mm  and  a  3.2  mm  diameter  collector  at 
air  speeds  of  3  and  5  m  s  As  shown  in  Apendix  III,  this 
increase  is  insufficient  to  have  any  effect  on  the  significance 
of  the  Jayaratne  and  Saunders  small  droplet  charging  experiments. 

The  results  presented  here  are  for  a  certain  set  of  experi¬ 
mental  parameters  sucn  as  collector  diameter,  cloud  droplet  spec¬ 
trum  and  impact  speeds.  The  net  contribution  of  the  rime  fea¬ 
thers  to  the  increase  in  collection  efficiency  must  be  obtained 
in  order  that  the  results  can  be  applied  to  a  different  experi¬ 
mental  situation.  Table  2  shows  the  proportion  of  the  droplets 
caught  to  the  total  number  in  the  air  flow  for  both  the  smooth 
ana  the  rough  collectors.  At  a  flow  speed  of  3  m  s  *  the  ratio 
increases  as  the  temperature  falls,  due  to  the  enhanced  growth  of 
the  rime  feathers  at  lower  temperatures.  Thus  the  volume  of 
droplets  collected  by  a  rimed  collector  may  be  increased  by  50% 
above  that  predicted  by  the  Lofflet  and  Muhr  equations  for  a 
smooth  collector.  However,  at  5  m  s~^,  the  ratio  for  the  rough 
collector  is  about  tne  same  as  that  of  the  smooth  collector. 
This  is  due  to  the  feather  growth  being  inhibited  at  higher 
impact  speeds  as  discussed  earlier. 

Col lect ion 

The  collection  efficiences  given  by  Loffler  and  Muhr  are 
only  valid  for  smooth  collectors.  When  tne  target  is  covered 
with  rime,  the  collection  efficiency  of  the  target  for  small 
droplets  does  increase,  particularly  at  low  temperatures.  It 
does  not,  however,  increase  with  the  velocity  because  feather 
growth  is  inhibited  wnen  the  rime  density  increases  with  increase 
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TABLE  1 


Collection  efficiency  (%)  given  by  Loffler  and  Muhr  for 
a  collector  diameter  of  5  and  3.2  mm  at  a  flow  speed  of 
3  m  s-1 

Droplet  5  mm  3.2  mm 
Diameter 


2.5 

0.00 

0.00 

5.0 

0.06 

0.21 

7.5 

0.66 

2.37 

10.0 

3.56 

10.94 

12.5 

11.28 

26.16 

15.0 

23.38 

40.94 

17.5 

35.84 

51.92 

20.0 

46.05 

59.81 

22.5 

53.87 

65.72 

25.0 

59.91 

70.36 

27.5 

64.73 

74.11 

30.0 

68.68 

77.21 

32.5 


74.81 


82.00 


TABLE  2 


Ratio  of  theoretical  volume  and  actual  volume 


to  total 

volume  in 

flow. 

Temp. 

3  m 

s-1 

5  i 

Smooth 

Rough 

Smooth 

-5 

0.300 

0.369 

- 

-8 

0.303 

0.430 

0.447 

-10 

0.304 

0.463 

0.448 

-12 

0.306 

0.468 

0.451 

-15 

0.308 

0.474 

0.453 

-20 

0.313 

0.482 

0.457 

collected 


Rough 

0.450 

0.452 

0.455 

0.460 

0.468 


in  tne  droplet  impact  velocity.  The  collection  efficiency  also 
increases,  up  to  a  certain  limit,  with  increasing  amounts  of  rime 
on  the  target.  The  collection  of  small  droplets  is  insufficient 
to  modify  the  soft-hail  stone  growth  rate,  and  hence  to  influence 
the  conclusions  of  Jayaratne  and  Saunders  (1985)  concerning  the 

sign  of  the  charge  transferred  when  ice  crystals  bounce  off  a 
soft-hailstone . 
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SECTION  5 


THE  COLLECTION  EFFICIENCY  OF  A  CYLINDRICAL  TARGET 
FOR  ICE  CRYSTALS 

Introduction 

The  collection  efficiency  of  an  ice  sphere  for  very  small 
ice  crystals  has  been  investigated  by  only  a  few  workers,  in¬ 
cluding  Hosier  and  Hallgren  (1960),  Hallgren  and  Hosier  (1960), 
Latham  and  Saunders  (1970)  and  Rogers  (1974).  There  have  been  no 
theoretical  studies,  owing  to  the  complexity  of  modelling  the 
movement  of  an  ice  crystal  through  a  fluid. 

Apart  from  Latham  and  Saunders,  all  the  previous  studies 
have  found  that  the  collection  efficiency  has  a  temperature 
dependence;  the  collection  efficiency  increases  as  the  tempera¬ 
ture  of  the  target  approaches  0°C  due  to  the  formation,  on  the 
surface,  of  a  pseudo-liquid  layer  which  enhances  the  adhesion 
efficiency.  Figure  17  summarises  their  results;  it  can  be  seen 
that  the  collection  efficiency  is  less  than  unity  in  all  cases. 
Rogers,  and  Latham  and  Saunders,  used  only  small  plate  crystals 
at  all  temperatures,  whereas  Hosier  and  Hallgren  dealt  with  the 
crystal  types  which  occur  at  the  temperature  at  which  the  experi¬ 
ments  were  performed.  Hosier  and  Hallgren  found  that  the  collec¬ 
tion  efficiency  was  a  maximum  of  about  0.2  at  about  -10°C,  de¬ 
creasing  to  0.05  at  -5°C.  They  suggested  that  this  change  in 
collection  efficiency  was  due  to  the  change  with  temperature  of 
ice  crystal  habit  from  plates  to  prisms. 

Latham  and  Saunders  (1970)  measured  the  collection 
efficiency  of  an  ice  sphere  of  diameter  2mm  for  plate  ice 
crystals  of  diameter  5  urn  at  an  impact  speed  of  3m  s-^.  The 
collection  efficiency  was  calculated  from  the  mass  of  ice 
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COLLISION  EFFICIENCY 


Figure  17 

Summary  of  previous  collision 
efficiency  work 


Latham  and  Saunders  (1970) 
Rodgers  (1974) 

Hallgren  and  Hosier  (1960) 
Hosier  and  Hallgren  (1960) 


crystals  collected  by  the  sphere  m  3  knowr,  tine.  The  collection 
efficiency  was  about  0.3  at  all  temperatures  1  r.  3  zero  electric 
field.  When  an  electric  field  of  about  1500  V  cm  ^  was  applied, 
the  collection  efficiency  increased  to  0.8.  They  state,  however, 
that  the  increase  is  due  to  an  increase  in  the  aggregation  effi¬ 
ciency,  and  not  in  the  collection  efficiency. 

Jayaratne  (1981)  found  that  the  continuous  formvar  replica¬ 
tor,  (Appendix  I),  has  a  collection  .fficiency  of  unity  for  ice 
crystals  greater  than  20pm  in  diameter  at  an  impact  speed  of 
4.7  m  s  ^ .  The  cross-sectional  shape  of  this  flow  was  a  rect¬ 
angular  jet,  and  so  the  collection  efficiency  for  a  20  pm  crystal 
may  not  be  unity  for  a  cylindrical  target.  Jayaratne  et  al . 
(1983)  calculated  the  Event  Probability  (the  collision  efficiency 
times  the  separation  probability)  by  comparing  the  charging 
current  to  two  targets  of  different  diameter  (see  Appendix  IV). 
The  event  probabilty  of  the  thinner  target  was  assumed  to  be 
unity  and  so  a  maximum  value  of  the  event  probability,  due  to  the 
collision  efficiency  not  being  known  at  that  time,  could  be 
calculated  for  the  thick  target.  The  charge  per  event  results 
were  not  corrected  for  the  collision  efficiency  of  the  target, 
and  so  will  tend  to  underestimate  the  charge  per  event  for  the 
smaller  ice  crystals. 

The  collection  efficiency  of  an  ice  crystal  for  water 
droplets  has  been  investigated  extensively  by  such  workers  as  Ono 
(1969),  Wilkins  and  Auer  (1970),  Harimaya  (1975)  and  Pitter 
(1977),  A  summary  of  this  work  is  shown  in  Figure  18.  The 
collection  efficiency  for  a  large  droplet  of  about  40  pm  diameter 
may  be  close  to  zero,  as  a  consequence  of  the  relative  velocity 
of  the  two  particles  being  close  to  zero.  This  section  describes 
the  work  performed  to  determine  the  collection  efficiency  of  a 
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cylindrical  target  for  ice  crystals  of  different  types.  Experi¬ 
ments  were  also  performed  to  determine  the  event  probability  of  a 
cylindrical  target. 


Experimental  Procedure 


The  method 

used 

was 

to 

compare  the  ice 

crystal 

number 

concentration  in 

the 

cloud 

to 

the  concentration 

found 

on  the 

cylindrical  target.  The  ratio  of  the  two  concentrations  is 
equivalent  to  the  target  collection  efficiency  and  will  increase 
with  increasing  flow  speeds  or  with  decreasing  target  diameter. 
The  experiments  were  performed  in  the  cloud  chamber  in  the  cold- 
room.  The  cloud  was  formed  by  nucleating  a  water  droplet  cloud 
in  the  usual  way.  The  size  and  number  concentration  of  the  ice 
crystals  formed  was  determined  from  formvar  covered  microscope 
slides  past  which  the  cloud  was  drawn. 

The  cylindrical  target  was  a  formvar  covered  glass  rod  of 
known  diameter.  These  coated  rods  were  prepared  in  the  same  way 
as  the  microscope  slides.  However,  before  exposure  to  the  crys¬ 
tal  cloud,  the  rod  was  dipped  in  cold  chloroform  to  activate  the 
formvar.  The  reason  for  this  is  that  the  ice  crystals  have  a 
finite  probability  of  bouncing  off  the  target  (see  Appendix  II), 
and  so  the  ice  crystal  concentration  would  be  underestimated. 
The  wet  formvar  ensures  that  the  target  has  a  sticking  efficiency 
of  unity,  and  that  all  ice  crystals  that  collide  with  the  target 
are  captured.  Thus  the  collection  efficiency  is  the  same  as  the 
collision  efficiency.  The  target  diameters  used  were  3,  4,  5  and 
6mm,  and  the  airflow  speeds  used  were  3,  4  and  5  m  s~*.  The 
airflow  speed  used  for  the  microscope  slide  was  7  m  s  so  the 
collection  efficiency  of  the  slide  was  assumed  to  be  unity  for 
all  but  the  smallest  (<  19  pm)  ice  crystals. 
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The  habit  of  ice  crystals  fall  into  three  distinct  types; 


columns,  plates  and  dendrites.  These  types  occur  at  around  -6,  - 
10  and  -15°C  respectively.  However,  small  dendrites  less  than 
80  pm  in  diameter  approximate  to  small  hexagonal  plates  as  the 
dendritic  arms  are  not  very  pronounced  at  this  size.  Thus  the 
experiments  were  performed  at  -6  and  -11°C,  with  ice  crystals 

corresponding  to  columns  and  hexagonal  plates.  The  ice  crystal 

4  7-3 

number  concentration  varied  from  10  to  10  m  and  the  ice 
crystal  diameter  ranged  from  10  pm  to  140  pm. 

Experiments  were  also  performed  with  wet  and  dry  formvar 
covered  glass  rods.  The  dry  target  simulates  a  hailstone  in  that 
the  sticking  efficiency  is  less  than  unity  and  ice  crystals  can 
collide  with  and  separate  from  the  target.  The  calculated  number 
concentration  was  again  compared  with  that  in  the  cloud,  to  give 
the  event  probability.  The  event  probability  could  be  compared 
with  that  previously  calculated  from  the  charging  currents  of  two 
simulated  hailstones  of  different  diameters,  as  described  in 
Appendix  IV. 

Measurement  of  the  Collision  Efficiency 

The  validity  of  the  assumpion  that  the  microscope  slide  has 
a  collection  efficiency  of  unity  was  confirmed  by  comparing  the 
number  concentration  on  the  5mm  rod  to  that  on  the  slide  for  an 
airflow  of  7ms  ^ .  Figure  19  shows  the  slide  number  concentra¬ 
tion,  Ns,  the  rod  number  concentration  N ,  and  the  ice  crystal 
diameter  . 

As  can  been  seen,  there  is  a  1:1  correlation  between  the  two 
concentrations  for  crystal  sizes  over  30  pm.  The  number 
concentration  on  a  dry  slide  and  a  wet  slide,  and  on  a  wet  rod 
and  dry  rod  were  also  compared.  Again,  there  was  a  1:1 
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relationship  between  t.ne  results. 

The  following  experiments  were  performed  with  wet  formvar 
targets  to  ensure  that  the  sticking  efficiency  was  unity. 

These  exeriments  were  performed  at  an  air  speed  of  3  m  s  ^ 
for  rod  diamters  of  3,  4,  5  and  6mm,  and  for  column  and  plate 
type  crystals  of  different  sizes.  Figures  20  and  21  show  the 
collision  efficiencies  against  ice  crystal  diameter  for  the 
different  rod  diameters,  for  plates  and  columns  respectively. 
The  collision  efficiency  of  the  3mm  rod  is  nearly  unity  for  all 
crystal  sizes.  It  can  be  seen  that  the  collision  efficiency  is 
unity  for  all  rods,  when  the  plate  diameter  is  greater  than  80  pm 
for  columns  of  length  100  pm. 

The  air  speed  was  varied  between  3  and  7  m  s  ^  for  the  5mm 
rod.  The  collision  efficiency  increases  with  increasing  air 
speed,  until  it  is  unity  for  all  crystal  diameters  examined  at 
7  m  s  S  Figures  22  and  23  are  for  plates  and  columns.  Again, 
the  collision  efficiency  for  a  column  type  crystal  is  lower  than 
that  of  a  plate  type  crystal  of  the  same  size. 

Measurement  of  the  Event  Probability 

In  these  experiments,  a  comparison  between  the  number  concen¬ 
tration  for  wet  and  dry  rods  was  made  in  order  to  obtain  the 
sticking  efficiency.  The  dry  rod  simulates  a  hailstone  whose 
sticking  efficiency  is  less  than  unity.  The  5mm  rod  was  used  so 
that  these  results  could  be  directly  compared  with  those  from  the 
charging  experiments.  Experiments  with  the  3mm  rod  were  also 
performed  to  investigate  the  effect  on  the  event  probability  of  a 
smaller  target  diameter.  The  event  probability  was  calculated 
from  the  product  of  collision  efficiency  and  one  minus  the 
sticking  efficiency. 
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F i u  r  e  2  f 


Collision  efficiency  of  Ci  Ff  irs"  t  r;  l3">»ter  rocs 
for  plate  crystals 


Xtal  DLam  (pm) 


F inure  21 


Collision  e‘','ici=-’ov  of  different 
diameter  rods  For  column  crystals 


Xtal  Dlam  (|jm) 


Figure  24  snows  tne  event  probability  against  ice  crystal 
diameter  for  Doth  rods  at-lO^C  an  air  speed  of  3ms  1 .  The  5mm 
target  has  a  lower  event  probability  that  the  3  mm  target  for  a 
particular  crystal  diameter.  The  event  probability  of  the  larger 
target  varies  from  0.1  to  0.35  for  the  crystal  diameter  range  10 

to  500  pm,  whereas  for  the  smaller  target  the  event  probability 

varies  from  0.25  to  0.55  over  the  same  range.  The  event  proba¬ 
bility  is  insensitive  to  crystal  type;  the  smaller  crystals  were 
columns,  plates  and  sectors,  whilst  the  larger  crystals  (>200  pm) 
were  sectors  and  dendrites.  The  event  probability  is  also  insen¬ 
sitive  to  cloud  and  target  temperature  but  is  very  sensitive  to 
impact  speed. 

Figure  25  shows  the  event  probability  for  both  targets 

against  crystal  diameter  for  an  air  speed  of  7  m  s-1.  The 
event  probability  of  the  5  mm  target  for  a  100  pm  crystal 

increases  from  0.1  at  3ms  1  to  0.3  at  7m  s_1.  Thus  there  is  a 
higher  proportion  of  the  ice  crystals  bouncing  off  the  taiget  at 
higher  impact  speeds. 

Discussion 

The  work  presented  here  shows  that  the  collision  efficiency 
of  a  cylindrical  target  for  ice  crystals  is  less  than  unity  if 
the  crystal  diameter  is  less  that  100pm.  Figures  19  tc  23  show 
that  the  collision  efficiency  is  dependent  on  crystal  diameter, 
D^,  speed,  v,  and  target  diameter,  D  .  These  paramet-rs  can  be 
related  by  the  following  equation: 

^  =  10vDx/Dc 

The  graph  of  ^  against  collision  efficiency  for  plates  and 
columns  is  shown  in  Figure  26a  and  b.  The  collision  efficiency 
of  a  target  for  plates  and  columns  cannot  be  related  directly  in 
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one  equation,  due  to  tne  different  deper.der,:^  on  crystal  shape. 
The  collision  efficiency  cannot  oe  related  t:  tae  equations  of 
Ranz  and  Wong  (1952)  or  those  of  Loffler  and  Muhr  (1977)  as  their 
equations  dealt  with  a  spherical  droplet,  and  the  results 
presented  here  show  that  a  hexagonal  plate  or  column  crystal 
cannot  be  idealised  to  a  sphere  as  assumed  by  Takahashi  (1978). 

The  above  relationship  was  used  to  calculate  the  collision 
efficiency  from  the  results  of  Latham  and  Saunders  (1970).  Their 
results  gave  a  collection  efficiency  of  0.3  and  the  equation 
gives  0.22,  which  is  in  fairly  good  agreement.  When  the  equation 
was  applied  to  the  results  of  Hosier  and  Hallgren,  it  gave  a 
collection  efficiency  of  0.5  as  opposed  to  0.1.  However,  sur¬ 
prisingly,  they  found  that  the  collection  efficiency  increased 
when  the  collector  diameter  increased  from  0.127mm  to  0.36mm,  and 
so  the  low  collection  efficiency  could  be  a  consequence  of  the 
very  small  collector  diameters  they  used. 

The  values  of  event  probability  calculated  here  are  in  good 
agreement  with  those  values  calculated  from  the  charging  experi¬ 
ments.  The  charging  determined  event  probability  against  collec¬ 
tion  efficiency  determined  event  probability  is  shown  in  Figure 
27.  These  results  presented  here  for  a  dry  target  will  tend  to 
overestimate  the  event  probability  slightly  at  warm  temperatures 
>  -10°C,  due  to  the  high  adheion  efficiency  of  a  riming 
ha ilstone . 

Concl us  ion 

The  collision  efficiency  of  a  target  for  ice  crystals  varies 
with  impact  speed,  crystal  diameter  and  rod  diameter.  The 
collision  efficiency  is  unity  for  all  crystal  sizes  greater  than 
100  pm.  The  results  presented  here  are  in  fairly  good  agreement 
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witn  those  of  previous  workers,  even  though  the  experimental 
procedure  was  radically  different.  The  event  probability  also 
varies  with  these  parameters,  with  the  probability  increasing 
with  increasing  impact  speed  and  crystal  diameter . 
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SECTION  6 


THE  CHARGING  OF  A  SIMULATED  SOFT-HAILSTONE  BY 
LARGE  ICE  CRYSTALS 


Introduction 

The  purpose  of  this  study  was  to  investigate  the  charging  of 
a  simulated  hailstone  under  a  range  of  cloud  conditions  with 
large  ice  crystals  greater  than  200  pm  in  diameter.  The  maximum 
size  achieved  by  Jayaratne  et  al.  (1983)  was  130  pm,  and  the 
present  work  is  designed  to  extend  the  charge  per  event-crystal 
diameter  dependence  and  the  velocity-crystal  diameter  dependence, 
determined  by  Jayaratne  et  al  .  to  larger  crystal  sizes.  Also  the 
charging  of  aircraft  was  simulated  by  impacting  ice  crystals  on  a 
rimed  and  an  unrimed  target  at  high  speeds  of  around  100  m  s-^. 

The  pioneering  work  of  Reynolds,  Brook  and  Gourley  (1957) 
proved  conclusively  that  a  soft  hailstone  can  become  electrified 
when  ice  crystals  bounce  off  it.  They  rotated  a  metal  sphere 
through  a  cloud  consisting  of  supercooled  water  droplets  and  ice 
crystals  at  -27°C,  and  measured  the  cnarge  acquired  by  the 
sphere.  They  detected  negative  charge  at  nigh  liquid  water 
contents,  and  positive  at  low.  Subsequent  work  by  Saunders  et 
al  .  (1984  >  showed  that  the  positive  charging  was  anomalous,  being 
due  to  the  high  rotational  forces  present  in  their  experimental 
set-up.  The  typical  liquid  water  content  was  atiout  2  g  m^,  and 
the  crystals  were  produced  m  concentrations  of  10 ;  m  ^  by 
seeding  the  cloud  with  dry  ice.  The  spheres  were  rotated  at  a 
high  speed  of  about  8  m  s  At  these  values,  they  measured  a 

charge  of  -160  fC  per  crystal  collision.  However,  their  esti¬ 
mates  of  crystal  concentration  may  have  been  in  error.  They  did 
not  detect  any  charging  when  the  sphere  was  rotated  in  water 
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dr opl e 's  3 1  one . 

Churcn  (  1966  !,  it,  a  similar  Set  of  nign  speed  experiments, 
found  tnat  a  naiistone  would  cnarge  positively  if  the  cloud 

temperature  was  warmer  than  -20 °C.  The  hailstone  charged 
negatively  if  the  vapour  supply  was  cut  off  during  a  run.  He 
measured  a  maximum  charge  per  event  of  1.6  fC. 

Takahashi  (1978)  also  repeated  the  Reynolds  experiments  at 
high  speed,  and  measured  negative  charge  at  high  liquid  water 
contents  below  -10°C  and  positive  at  low  values  of  liquid  water 
content.  The  charge  per  collision  ranged  from  -16  fC  to  +30  fC 
for  a  100  pm  diameter  crystal . 

Jayaratne  et  al.  (1983)  provided  the  first  comprehensive 
investigation  into  the  charging  of  a  soft  hailstone  by  the  coll¬ 
ision  of  ice  crystals  at  velocities  appropriate  to  those  found  in 
thunderstorms.  When  they  rotated  a  target  through  a  cloud  of 
supercooled  water  droplets,  riming  occurred  on  the  leading  edge 
but  no  charging  to  the  hailstone  was  detected.  The  cloud  was 
then  seeded  with  a  wire  wmch  had  been  cooled  to  liquid  nitrogen 
temperature,  and  ice  crystals  grew  in  the  supersaturated  environ¬ 
ment.  Strong  positive  cnarging  was  detected  at  a  temperature  of 
about  -18°C  for  liquid  water  contents  greater  tnan  0.3  a  m 3  ; 
negative  charging  was  detected  below  this  'reversal'  temperature  , 
typically  -180C,  unless  the  liquid  water  content  was  unrealis¬ 
tically  high.  Tne  charge  per  event  for  a  125  u~  diameter  crystal 
was  about  +10fC  at  an  impact  speed  of  3  m  »  ^  at  a  temperature  of 
—  1 0 ° C .  Tne  present  work  extends  these  results  to  larger  crystals 
and  higher  velocities. 
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Experimental  Procedure 

Tne  low  speed  exer  iments  were  performed  entirely  witnm  t.ne 
cloud  cnamber  in  tne  cold-room.  The  targets  were  mounted  on  the 
rotating  apparatus,  and  could  be  moved  through  the  cloud  with 
their  axes  normal  to  the  direction  of  motion.  The  target  speed 

was  typically  3  m  s  and  riming  occurred  on  the  leading  edge 

The  target  was  connected  to  an  electrometer  which  had  a  sensi¬ 
tivity  of  1  mV  =  10  ^  A  which  could  detect  the  charge  to  the 
target  if  the  charging  rate  was  greater  than  0.2  pC  s-^.  The 
output  was  recorded  on  a  chart  recorder.  Two  targets  were  used; 

one  of  diameter  5mm  and  the  other  of  diameter  0.5mm.  They  were 

placed  at  opposite  arms  of  the  apparatus,  as  shown  in  Figure  2. 
The  event  probability  (see  Appendix  IV  and  Section  5)  is  assumed 
to  be  unity  for  the  smaller  target. 

The  cloud  chamber  consisted  of  a  metal  box  divided  into  two 
parts;  the  upper  was  the  initial  growth  chamber  for  the  ice 
crystals,  and  the  lower  was  the  experimental  chamber  where  the 
ice  crystals  interacted  with  the  simulated  hailstone.  A 
supersaturated  environment  was  maintained  in  Doth  chambers;  The 
cloud  m  the  upper  chamber  was  seeded,  and  the  crystals  grew  to 
about  200  pm.  At  this  point,  the  trapdoor  was  opened  and  the 
crystals  fell  through  to  the  bottom  chamber.  After  about  3 
minutes,  the  crystals  were  about  600  pm  in  diameter.  The  maximum 
growth  rate  occurred  at  -15°C  with  dendritic  crystals;  the 
maximum  sized  plate  crystal  was  about  300  pm.  There  was  an 
unavoidable  temperature  gradient  of  1°C  m~*  over  the  height  of 
the  cloud  chamber  but  the  region  where  the  riming  occurred  was 
well  mixed  and  so  had  uniform  temperature. 

The  ice  crystal  number  concentration  was  measured  by  drawing 
a  sample  of  cloud  past  a  dry  formvar  covered  microscope  slide  at 
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7  m  s'*.  The  slide  was  then  sprayed  with  cold  chloroform  and 
allowed  to  dry  before  being  analysed  under  a  microscope.  The 
actual  concentration  was  calculated  from  the  equations  described 
in  Appendix  I.  The  collision  efficiency  of  the  target  for  small 
ice  crystals  was  taken  into  account,  using  the  results  presented 
in  Section  5.  The  Rime  Accretion  Rate  was  determined  from  the 
elevation  of  the  rime  surface  temperature  above  ambient,  caused 
by  the  release  of  heat  when  the  droplets  collide  with  and  freeze 
on  the  target,  as  detailed  in  Appendix  II. 

Thus,  for  each  experimental  run,  the  charging  current  to 
both  simulated  soft-hailstones,  the  rime  accretion  rate,  the 
cloud  temperature  and  the  size  and  number  concentration  of  the 
ice  crystals  could  be  measured. 


seeded  between  1  and  2  minutes  previously.  Jayaratne  et  al . 
noted  an  initial  positive  peak  when  crystals  first  hit  the  riming 
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•.  i  r  ^  ■.  ,  w~.  i  :n  occur  red  irrespective  or  tne  cloud  temperature  . 

1ms  whs  ice  t  u  a  Mgn  concentration  of  small  ice  crystals  soon 
after  seeding  trie  cloud.  In  tne  present  case,  tnere  was  only  a 
uniform  cloud  of  larger  crystals  and  so  the  positive  peak  was  not 
detected.  The  positive  peak  could  be  reproduced  if  the  cloud  was 
seeded  at  the  bottom  of  the  chamber  near  the  target  apparatus. 
The  shape  of  the  charging  current  and  the  variation  of  the  liquid 
water  content  with  time  was  very  similar  to  that  of  Jayaratne. 
The  maximum  diameter  dendrite  produced  in  the  chamber  was  800  pm, 
and  the  charge  per  event  for  this  sue  of  crystal  was  +220  fC. 

The  experiments  were  repeated  at  -6°C  where  the  predominant 
ice  crystal  habit  is  columns  and  prisms.  Tne  time  surface  tem¬ 
perature  was  around-5°C.  The  maximum  reproducible  column  length 
was  380  pm  which  gave  a  charge  per  event  of  100  fC.  Figure  29 
shows  typical  charging  currents  to  both  targets  at  this  tempera¬ 
ture.  The  event  probability  for  this  run  witn  300  pm  columns  was 
calculated  to  be  0.25,  which  is  a  slightly  smaller  value1  than  at 
lower  temperat  jr-s  di-.  to  the  nign~r  idr.esion  ef  f  ici-ncy  at 
warmer  tempe r a t u r e s .  figure  30  snows  tne  variation  with  time  of 
the  ice  crystal  size  and  number  c.T.cV.nt  ration,  and  the  rime 
•accretion  rate.  The  largest  plat®  crystal  obtainable  was  300  pm,, 
leading  to  a  charge  per  event  of  80  f C. 

Figure  31  snows  tne  charge  per  event  against  crystal  dia¬ 
meter.  Jayaratne  et  al.  found  a  relationship  between  chat  >'  and 

4 

size  oi  the  form  q  «.  d  for  crystals  up  to  about  120  pm;  in  -  i  ? :<_■ . 

4 

Their  result  was  confirmed  for  small  crystals  a  ,d  the  d  line  is 
shown  in  the  figure.  For  crystals  greater  than  150  pm  diameter, 
the  charge  per  event  is  less  than  that  predicted  by  extrapolating 
the  d*  relationship  to  larger  crystal  sizes.  Above  450  pm  there 
is  only  a  slight  increase  in  the  charge  transferred  with  an 
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The  Negative  Charging  Regime 

At  temperatures  below  the  reversal  temperature  of  -16°C,  the 
charging  to  the  target  is  predominantly  negative.  The  main 
crystal  habit  is  plates  and  bullets  at  these  temperatures. 
Experiments  were  performed  to  measure  the  negative  charge  per 
event  for  large  ice  crystals. 


Figure  32  snows  typical  runs  at  -25°C  and  -29°C;  the 


charging  is  negative  throughout  the  run.  The  predominant  crystal 
type  was  bullet  clusters,  which  separated  more  charge  than  a 
single  crystal  of  the  same  size.  This  is  reflected  in  the  magni¬ 
tude  of  the  charging  current.  When  they  first  appear  at  the 
target,  the  clusters  are  about  400  pm  in  diameter  and  reach  about 


1000  pm  in  diameter  at  the  end  of  the  run. 


The  results  of  Jayaratne  et  al .  were  confirmed  by  seeding 
the  cloud  in  the  lower  chamber.  The  initial  positive  peak,  was 
observed,  followed  by  negative  charging  during  the  rest  of  the 
run.  Figure  33  shows  such  a  run  at  -21°C;  the  maximum  crystal 
diameter  was  140  pm.  The  figure  also  shows  the  crystal  size  and 
numo? r  concentration  during  the  run.  No  bullet  clusters  were 
observed,  only  plates. 


va . *s  o£  charge  per  event  from  these  experiments  were 
to  those  found  by  Jayaratne  for  small  crystals.  Figure 
~na  rge  per  event  against  crystal  diameter;  the  line 
r:  : rom  the  Jayaratne  line  at  a  single  crystal  diameter 
Tie  bullet  clusters  separated  more  charge,  as 
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shown  in  the  figure.  A  possible  explanation  for  this  is  dis¬ 
cussed  later  . 

The  Effect  of  the  Rime  Accretion  Rate  on  the  Charge  per  Event 

A  parameter  which  controls  the  sign  and  the  magnitude  of  the 
charge  per  event  is  the  rime  accretion  rate.  At  low  values,  the 
charging  to  the  target  is  negative  at  all  temperatures,  whereas 
at  high  values  the  charging  is  positive.  Jayaratne  et  al .  found 
that  the  sign  and  magnitude  of  the  charge  per  event  was  dependent 
on  the  cloud  liquid  water  content,  which  can  be  calculated  from 
the  rime  accretion  rate  if  the  droplet  spectrum  and  the  collec¬ 
tion  efficiencies  of  the  target  for  the  droplets  is  known.  The 
method  employed  in  this  study  was  to  measure  the  rime  accretion 
rate,  directly,  from  the  temperature  elevation  of  the  target 
above  ambient,  which  is  caused  by  the  freezing  of  the  water 
droplets  (see  Appendix  II).  Thus  the  collection  efficiencies  and 
the  droplet  spectrum  need  not  be  known. 

In  the  experiments  presented  here,  the  boiler  power  was 

varied  so  that  the  rime  accretion  rate  varied  from  0  to 
-2  -1 

25  mg  cm  min  at  different  temperatures  m  different  runs. 

The  graph  shown  in  Figure  35  illustrates  the  dependence  of  the 

sign  of  the  charge. per  event  on  the  rime  accretion  rate  versus 

temperature.  The  dependence  of  the  magnitude  of  the  charge  per 

event  for  a  150  pm  crystal  on  the  rime  accretion  rate  at  -10°C  is 

-  2  -1 

shown  in  Figure  36.  At  3  mg  cm  min  the  charge  per  event  is  - 

-2  -1 

2fC  rising  to  70  fC  at  a  rime  accretion  rate  of  25  mg  cm  min 


The  Effect  of  Ice  Crystal  Growth  Temperature 
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liquid  water.  The  maximum  growth  rate  occurs  with  dendritic 
crystals  at  -15°C  and  so  the  top  chamber  was  warmed  to  -15SC 
whilst  keeping  the  lower  chamber  well  below  the  reversal  tempera¬ 
ture.  Surprisingly,  the  resulting  charge  transfer  to  the  target 
was  positive,  even  though  the  lower  chamber  temperature  was  below 
-20°C.  The  sign  of  the  charging  is  opposite  to  that  normally 
obtained  at  these  temperatures.  Figure  37  shows  two  runs;  one 
at  -28°C  with  the  upper  chamber  at  -29°C  and  the  other  at  -28°C, 
with  the  upper  chamber  temperature  of  -15° C.  The  two  charging 
currents  are  of  opposite  sign,  even  though  the  air  and  rime 
surface  temperatures  of  the  two  runs  are  the  same.  The  rime 
accretion  rate  for  both  runs  is  also  similar. 

To  investigate  this  unusual  behaviour  for  crystals  grown  at 
-15°C  but  interacted  at  colder  temperatures  where  negative  target 
charging  is  expected,  these  experiments  were  repeated  using 
stationary  targets.  Figure  38  shows  the  charging  currents  to  the 
targets  when  the  airflow  is  turned  on,  then  off.  The  top  chamber 
temperature  was  -15 °C,  whilst  the  lower  chamber  temperature  was 
-28°C.  The  pump  control ing  the  crystal  flow  past  the  target  was 
turned  on  at  certain  times  to  measure  the  charging  current  in 
order  to  prevent  the  tube  containing  the  target  from  riming 
closed.  At  time  t=0,  the  trapdoor  was  opened  and  the  pump  was 
activated  at  t=20s;  the  sign  of  the  charging  current  was  posi¬ 
tive  even  though  the  interactions  were  occurring  below  the  rever¬ 
sal  temperature.  The  kettle  was  then  switched  off  and  the  crys¬ 
tals  depleted  the  cloud  and,  at  t=90s,  the  sign  was  negative.  At 
this  point,  the  crystals  were  evaporating  in  the  unsaturated 
environment.  Vapor  was  reintroduced  at  t=240s;  shortly  after¬ 
wards  the  pump  was  turned  on  and  the  sign  was  again  positive;  at 
this  point  the  liquid  water  content  was  estimated  to  be  about 


0.5  g  .  The  above  procedure  was  repeated  witn  upper  and  lower 
chamber  temperatures  of  -23  and  -28°C  respectively,  wnen  the 
target  charged  negatively  (Figure  39).  The  small  positive  peaic 
at  the  end  of  the  run  is  due  to  the  very  high  liquid  water 
content  of  3-4  g  m  ^  present  in  the  chamber. 

The  initial  positive  charging  current  was  reproducible  if 
the  upper  chamber  temperature  was  above  the  reversal  temperature; 
below  -18°C  the  charging  was  always  negative.  Negative  charging 
for  a  warm  upper  chamber  could  also  be  obtained  if  the  liquid 
water  content  was  below  0.1  g  m-^.  The  consequences  of  these 
results  are  discussed  later. 
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Chargim  in  the  P resence _ of _ Ice  Crystals  Alone 

A  cloud  consistir.q  of  ice  crystals  only  *'i.;  retained  by 
seeding  a  supercooled  cloud  of  water  droplets,  switching  off  tne 
vapor  supply  and  allowing  tne  crystals  to  use  up  all  the  availa¬ 
ble  water  vapor.  In  this  manner,  experiments  could  be  performed 
without  any  water  droplets  being  present.  The  apparatus  was 
rotated  at  3  m  s  *  through  a  cloud  of  120  pm  plate  crystals;  no 
charging  was  detected  above  the  noise  on  the  amplifier  which  was 
consistent  with  the  results  of  Jayaratne  et  al.  However,  when 
ice  crystals  greater  than  200  pm  in  diameter  interacted  at  a 
speed  of  3  m  s  ^  negative  charging  was  measured.  Figure  40  shows 
the  charge  per  event  against  crystal  diameter.  The  magnitude  of 
the  charge  per  event  ranges  from  -0.1  fC  for  a  200  pm  crystal  to 
-1  fC  for  a  500  pm  crystal.  These  values  are  about  1  order  of 
magnitude  smaller  than  the  values  of  negative  charging  obtained 
at  low  liquid  water  content,  and  two  orders  of  magnitude  lower 

than  the  positive  charging  obtained  at  -15°C  at  a  liquid  water 

,  ,  -3 

content  of  1  g  m 

Low  Velocity  Charging 

Jayaratne  et  al .  (1983)  found  that  tne  charge  per  event 

increased  with  impact  velocity  raised  to  the  power  3.2  for  small 
ice  crystals.  Similar  experiments  were  performed  with  200  pm  ice 
crystals  in  a  mixed  cloud  at  -12°C  and  -25°C.  Tne  maximum  per¬ 
missible  velocity  with  the  rotating  apparatus  was  6  m  s  The 

cloud  liquid  water  content  was  controlled  to  give  the  same  rime 
accretion  rates  for  all  the  velocities  investigated.  Figures  41 
and  42  show  the  charge  per  event  at  different  velocities  at  the 
two  temperatures.  The  figures  also  show  the  results  discussed  in 
the  following  section  at  10  and  25  ms"*  . 
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High  Velocity  Charging 


Charging  with  Ice  Crystals  Alone 

In  these  experiments,  the  amplifier  ratio  was  changed  so 
12 

that  lmV  S  10  A.  At  the  greater  sensitivity  the  amplifier  was 

swamped  by  the  large  charges  being  separated.  The  crystal  cloud 

was  produced  in  the  usual  way,  and  the  crystals  were  allowed  to 

4  -3 

fall  out  until  only  a  small  concentration  of  about  10  m 
remained. 

The  pure  crystal  cloud  was  drawn  past  a  5mm  diameter  metal 
target  at  a  speed  of  110  m  s  The  crystal  concentration  was 

determined  from  a  microscope  slide  just  before  the  pump  was 
switched  on.  It  was  assumed  that  the  crystal  concentration  drawn 
past  the  target  for  the  first  15  seconds  of  a  run  was  constant, 
after  which  time  the  crystal  concentration  was  depleted.  The 
charge  per  event  was  calculated  frcm  the  crystal  concentration 
and  the  initial  negative  peak.  The  event  probability  was  found 
to  be  unity,  as  there  was  no  evidence  of  any  crystals  on  the 
target  at  the  end  of  a  run.  Strong  negative  charging  was  mea¬ 
sured  at  all  temperatures  and  all  crystal  diameters. 

Figure  43  shows  the  charge  per  event  against  crystal  dia¬ 
meter  at  a  speed  of  110  m  s-1.  In  agreement  with  Jayaratne  et 
al . ,  the  sign  of  the  charging  current  is  negative  when  the  target 
is  evaporating.  The  figure  illustrates  that  large  amounts  of 
charge  are  being  separated?  for  example,  a  500  crystal  sep¬ 
arates  about  -10  pC,  which  is  3  orders  of  magnitude  greater  than 
the  charge  measured  at  3  m  s-^. 

These  experiments  were  repeated  with  a  previously  rimed 
target.  The  target  was  rimed  for  a  few  seconds;  the  cloud  was 
then  seeded  and  the  experiments  conducted  as  beFore  at  110  m  s'1. 
There  was  no  effect  due  to  the  presence  of  the  rime,  as  illus- 
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trated  by  Figure  44  wmch  plots  points  from  the  two  sets  of 
experiments  at  equal  values  of  crystal  sizes.  There  is  a  close 
1:1  relationship  between  the  charge  per  event  of  the  rimed  tar¬ 
get,  QR,  with  that  of  the  unrimed  target,  Q. 

Charging  in  a  Mixed  Cloud 

The  riming  experiments  were  conducted  with  the  stationary 
target  by  allowing  a  few  droplets  to  remain  in  the  ice  crystal 
cloud.  The  Macklin  and  Payne  equation  (Appendix  II)  shows  that 
at  these  high  velocities,  a  very  low  liquid  water  content  of 
0.1  g  m  ^  will  lead  to  a  temperature  rise  from  which  a  rime 
accretion  rate  of  14  mg  ^  min  ^  may  be  calculated. 

The  apparatus  was  modified  so  that  the  cloud  could  simul¬ 
taneously  be  drawn  past  three  targets  at  different  speeds,  which 
enabled  direct  comparison  between  three  speeds  under  the  same 
cloud  conditions  to  be  made.  The  target  diameter  was  1.6  mm  and 
the  flow  speeds  used  initially  were  75,  50  and  25  m  s  ^ .  The 
event  probability  is  unity  for  all  crystal  diameters  at  these 
velocities.  Figure  45  shows  currents  from  the  three  targets  at  a 
temperature  of  -15°C.  The  fact  that  the  shape  of  the  curves  is 
the  same  verifies  that  the  targets  are  experiencing  the  same 
cloud  conditions.  Figure  46shows  positive  charge  per  event 
against  crystal  diameter  for  these  three  velocities.  It  can  be 
seen  that  the  three  lines  are  roughly  parallel  but  that  more 
charge  is  separated  at  the  lower  speed. 

The  apparatus  was  also  used  at  a  temperature  of  -25°C.  The 
sign  of  the  charging  was  negative  but  the  same  pattern  of  velo¬ 
city  dependence  was  observed.  Figure  47  shows  the  charging 
current  to  the  targets  at  this  temperature.  The  shape  of  the 
current  versus  time  is  the  same  for  all  the  targets,  thus  con- 
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Figure  UU 

q  vs  QR  at  110  m/s 


12 


QR  (fC) 


firming  that  they  are  experiencing  the  same  cloud  conditions. 
Figure  48  shows  the  charge  per  event  against  crystal  diameter, 
with  maximum  values  at  25  m  s~^.  The  gradients  of  the  graphs  for 
both  positive  and  negative  charging  is  0.7.  A  possible  explana¬ 
tion  for  the  velocity  dependence  is  discussed  in  the  following 
section. 

Discussion 

Large  amounts  of  charge  are  generated  when  an  ice  crystal 
hits  and  rebounds  from  a  riming  hailstone.  The  magnitude  and 
sign  of  the  charge  is  dependent  on  the  rime  accretion  rate,  the 
ice  crystal  diameter,  the  impact  speed  and  the  cloud  temperature. 
Jayaratne  et  al .  (1983)  found  the  dependence  of  the  charge  per 
event  with  crystal  diameter  up  to  120  pm.  These  positive 
charging  results  were  extended  to  a  diameter  of  800  pm  and  both 
sets  of  results,  whicn  have  been  corrected  with  revised  values  of 
event  probability  determined  in  the  present  study,  are  shown  in 
Figure  49.  It  can  be  seen  that  the  charge  per  event  changes 
shape  at  150  pm  and  480  urn.  For  example,  a  350  pm  crystal, 
according  to  the  Jayaratne  prediction,  would  separate  +1.25  pC, 
whereas  these  experiments  have  shown  that  it  only  separates  +0.1 
pC,  roughly  one  order  of  magnitude  less.  This  would  suggest  that 
there  is  some  mechanism  which  is  limiting  the  amount  of  charge 
being  separated.  There  was  no  dependence  of  the  charge  per  event 
on  crystal  type. 

In  the  negative  charging  regime,  the  deviation  from  the 
Jayaratne  experimental  line  occurs  at  210  pm;  Figure  50  shows 
both  the  new  and  the  Jayaratne  data.  The  deviation  from  the 
previous  data  is  greater  than  for  the  positive  case;  for  example, 
a  350  pm  crystal,  from  the  Jayaratne  data,  is  predicted  to  sep- 
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3  m/ s  with  revised  Jayaratne  data  for  the 
negative  charging  regine 


□  :  clusters  at  -25DC 


The 


arate  57  t'C  but,  in  reality,  it  only  separates  1.8  fC. 
largest  crystal  produced  in  this  regime  was  500  um.  The  charge 
per  event  for  the  bullet  clusters  was  greater  than  that  for  a 
single  crystal  of  the  same  size  due  to  the  cluster  shattering  on 
impact . 

The  high  speed  collision  experiments  were  designed  to  simu¬ 
late  the  charging  of  an  aircraft  as  it  flies  through  a  cloud.  In 
regions  of  zero  liquid  water  content,  the  target  charged  nega¬ 
tively  at  all  temperatures.  However,  the  collection  efficiency 
for  water  droplets  at  these  speeds  is  such  that  a  very  low  liquid 
water  can  give  a  high  rime  accretion  rate  and  consequently  lead 
to  positive  charging.  This  was  the  case  with  the  1.6mm  targets; 
under  cloud  conditions  which  gave  negative  charging  to  the  5mm 
target,  the  smaller  target  measured  positive  charging  due  to  its 
greater  collision  efficiency.  Thus,  it  was  quite  difficult  to 
ensure  a  zero  liquid  water  content  in  the  crystal  only  experi¬ 
ments  with  the  1.6mm  targets.  The  velocity  dependence  was  oppo¬ 
site  to  that  expected,  in  that  the  high  velocities  produced  less 
charge  than  at  lower  velocities.  Figure  51  illustrates  the 
charge  per  event  for  a  200  pm  crystal  against  air  speed.  It 
shows  that  the  charge  per  event  is  a  maximum  at  about  25  m  s  . 
This  can  be  explained  in  terms  of  the  contact  time  and  the 
sticking  efficiency.  As  the  contact  time  between  the  hailstone 
and  the  ice  crystal  increases,  the  charge  trarsfer  also  increases 
until  some  equilibrium  between  the  particles  is  reached  and  the 
charge  transfer  decreases.  Thus,  if  the  contact  time  for  an 
impact  speed  of  110  m  s”^  is  less  than  that  at  25  m  s  then 

less  charge  would  be  separated  at  the  higher  speed. 
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SECTION  7 


DISCUSSION,  CONCLUSION  AND  FUTURE  WORK 

The  hot-wire  device  discussed  in  Section  3  is  useful  for 
measurements  of  the  liquid  water  content  in  clouds.  The  novel 
approach  here  is  that  the  wire  is  heated  by  a  chain  of  fixed 
amplitude  voltage  pulses  whose  mark-space  ratio  is  varied  to 
control  the  mean  power  dissipation  in  the  wire.  A  wire 
temperature  of  76°C  was  chosen  for  maximum  sensitivity  and  to 
avoid  the  formation  of  an  insulating  layer  of  liquid  over  the 
wire  surface  at  high  temperatures.  The  output  is  simply  cali¬ 
brated  using  a  cloud  of  known  liquid  water  content  (determined  by 
another  longer  and  more  tedious  method).  The  device  is  accurate 
to  4%  at  1  g  m-3. 

A  cylindrical  target  has  a  collection  efficiency  for  water 
droplets  which  is  dependent  on  the  target  diameter,  the  air 
temperature  and  the  droplet  diameter.  Theoretical  collision 
efficiencies  have  been  calculated  for  smooth  targets;  however, 
rough  targets  present  a  more  difficult  problem.  The  experiments 
described  in  Section  4  show  that  the  collision  efficiency  is 
increased  at  temperatures  below  -10°C  and  at  air  speeds  below  5 
ms  3  due  to  the  growth  of  'rime  feathers'  on  the  target  surface. 
At  higher  speeds,  the  growth  of  the  feathers  is  inhibited.  An 
equivalent  collector  diameter  has  been  defined  as  the  diameter  of 
a  smooth  collector  whose  collection  efficiency  is  the  same  as 
that  of  the  rimed  collector.  At  temperatures  below  -10°C  at  an 
air  speed  of  3  m  s  3  the  equivalent  collector  diameter  of  a  5mm 
rimed  target  is  3.2mm.  At  5  m  s-3  for  the  same  conditions  the 
diameter  is  4.8mm.  At  the  lower  speed,  the  weighted  collection 
efficiency  is  increased  by  almost  50%  due  to  the  presence  of  the 
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rime  feathers. 


The  collision  efficiency  of  a  cylindrical  target  for  ice 
crystals  has  not  been  modelled  theoretically  and  has  only  been 
studied  in  a  limited  way  before.  The  results  described  in  Sec¬ 
tion  5  show  that  the  collision  efficiency  is  unity  for  crystals 
greater  than  80  um  at  a  flow  speed  of  3  m  s~*.  A  simple  rela¬ 
tionship  between  crystal  diameter,  target  diameter  and  flow  speed 
was  found. 

The  target,  moving  at  3m  s  ^ ,  did  not  undergo  detectable  or 
significant  charge  transfer  in  a  cloud  consisting  only  of  crys¬ 
tals  if  they  were  smaller  than  120  pm.  Above  this  size,  a  weak 
negative  charge  was  measured  which  was  dependent  on  crystal  size. 
At  high  speeds  strong  negative  charging  of  about  1  pC  per  event 
was  measured  at  all  temperatures.  The  charge  per  event  was  found 
to  be  dependent  on  the  impact  speed  with  a  maximum  value  at 
around  25  m  s  Impacts  on  airplanes  at  high  speeds  may,  there¬ 
fore,  not  lead  to  charge  transfer. 

The  above  mentioned  work  was  all  preliminary  to  the  achieve¬ 
ment  of  the  main  objectives  of  this  study.  In  order  to  calculate 
the  charge  separated  by  an  individual  crystal  bouncing  off  an  ice 
target  in  the  laboratory  experiments,  the  concentration  of  crys¬ 
tals  in  the  cloud  is  needed,  the  collection  efficiency  or,  in  our 
case,  the  event  probability  is  required,  and  the  current  to  the 
ice  target  must  be  measured  to  provide  a  value  of  the  charge  per 
event . 

Experiments  carried  out  have  followed  the  lines  of  previous 
investigations,  with  the  exception  that  larger  ice  crystals  than 
have  been  available  before  have  been  used.  Studies  of  charge 
transfer  for  both  the  positive  and  negative  regimes  at  tempera¬ 
tures  on  either  side  of  the  reversal  temperature,  have  been  made 
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with  ice  crystals  up  to  800  pm  in  diameter.  Previous  results 
with  crystals  up  to  100  pm  had  shown  that  the  charge/size  depen¬ 
dence  was  to  the  fourth  power.  This  result  has  been  verified. 
However,  at  larger  sizes  the  charge  transfer  power  dependence 
decreases  dramatically,  falling  to  a  value  of  order  unity  for 
crystals  around  500  pm.  This  decrease  with  larger  crystal  sizes 
is  noted  for  both  positive  charging  and  negative  charging,  which 
suggests  that  the  charge  transfer  mechanism  is  the  same  for  both 
signs  of  charging.  With  crystals  alone  in  the  cloud  the  charge 
transfer  was  weak  and  this  situation  is  probably  not  of  signifi¬ 
cance  in  thunderstorms. 

4 

Jayaratne  and  Saunders  (1985)  extrapolated  their  d  rela¬ 
tionship  to  larger  sizes  than  those  investigated  in  the  labora¬ 
tory  experiments  to  account  for  lower  positive  charge  centers 
(LPCC).  They  calculated  that  a  3mm  diameter  graupel  particle 
falling  1  km  through  a  crystal  concentration  of  104  m~3  would 
acquire  a  charge  of  +35  pc  if  the  charge  separated  in  each  in¬ 
teraction  was  500  fc.  They  assumed  that  the  size  of  the  crystals 
was  200-300  pm  with  an  event  probability  of  0.1.  However,  the 
results  presented  here  show  that  a  300  pm  crystal  will  separate 
aciout  80  fC  and  that  a  3mm  graupel  particle  has  an  event  proba¬ 
bility  of  0.5  at  3ms  .  Thus  the  particle  will  acquire  about 
29  pC.  Surprisingly,  therefore,  the  lower  than  expected  values 
of  charge  transfer  for  larger  crystals  have  not  led  to  a  signifi¬ 
cant  fall  in  the  total  charge  transfer,  simply  because  of  the 
increased  values  of  event  probability  with  the  larger  crystals. 

The  charge  separated  when  an  ice  crystal  of  diameter  d  re¬ 
bounds  from  a  hailstone  can  be  calculated  from  the  following 
relationship:  Q  =  Adav^  where  v  is  the  impact  velocity. 
Jayaratne  et  al.  found  that  the  constants  a  and  b  take  the  values 
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4  and  3  respectively.  The  present  study  has  shown  that  b=2  for 
low  velocities  and  a  takes  values  to  between  0.3  and  2.  The 
range  of  values  of  a  is  due  to  the  dependence  on  crystal  size. 
Calculations  are  now  under  way  which  include  these  latest  results 
in  a  thunderstorm  model,  so  that  the  rate  of  magnitude  of  elec¬ 
trification  may  be  compared  with  the  field  data. 

A  new  theory  of  the  charge  transfer  mechanism  is  being 
considered.  A  vast  array  of  experiments  have  now  been  performed 
using  a  range  of  experimental  conditions  and  any  theory  must  be 
consistent  with  the  results.  Recently,  Baker  et  al .  (1987) 

suggested  that  the  charge  transfer  is  associated  with  the  rela¬ 
tive  growth  rates  of  the  interacting  particles;  effectively,  the 
surface  which  grows  from  the  vapor  fastest,  charges  positively. 
However,  the  present  results  suggest  that  it  is  the  history  of 
the  growth  process  that  is  important  -  for  example,  crystals 
grown  at  -15®C  and  interacted  at  -28°C  charged  the  target  posi¬ 
tively;  whereas,  at  -28°C,  negative  charging  is  expected.  A 
possible  explanation  that  is  being  investigated,  is  that  the 
freezing  time  of  droplets  on  the  target  controls  the  concentra¬ 
tion  of  dislocations  in  the  ice  lattice  and,  similarly,  the  ice 
crystals  include  more  dislocations  in  their  surface  when  they 
grow  rapidly  from  the  vapor.  A  difference  in  dislocation  concen¬ 
trations  on  the  two  interacting  surfaces  leads  to  charge  transfer 
because  the  dislocations  are  charged.  Evidence  for  charged  dis¬ 
locations  is  available,  but  more  work  is  needed  before  this 
concept  can  be  assumed  to  be  the  charge  transfer  mechanism. 

The  size  dependence  results  need  to  be  extended  to  still 
larger  crystal  sizes,  preferably  up  to  l-2mm  diameter  crystals, 
to  discover  if  the  observed  dependence  continues.  Also,  interac¬ 
tions  between  millimetre-sized  graupel  particles  having  had 
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different  growth  histories  may  be  important.  One  method  of 
growing  graupel  in  the  cloud  chamber  would  be  to  introduce  some 
form  of  recirculation  into  the  cold  room  in  order  to  extend  the 
crystal  lifetime  and  to  promote  riming  of  the  ice  crystals. 

The  moving  target  experiments  presented  here  simulate  a 
soft-hail  stone  moving  through  a  cloud;  however,  the  height  of 
the  cold-room  would  allow  ice  spheres  of  different  diameters  to 
be  dropped  through  a  mixed  cloud  in  a  realistic  manner.  A  series 
of  induction  rings  may  be  arranged  vertically  to  measure  the 
charge  on  the  riming  hailstone  as  it  falls.  It  can  be  calculated 
that  the  hailstone  will  be  falling  at  its  terminal  velocity  in 
the  bottom  half  of  the  cloud  chamber,  and  given  typical  crystal 
concentrations  and  sizes  the  charge  acquired  by  the  hailstone 
will  be  in  the  region  of  1  pC.  A  further  extension  of  these 
experiments  would  be  to  introduce  an  electric  field  into  the 
chamber,  thus  allowing  inductive  ice-ice  collisions  and  inductive 
water-ice  collisions  to  be  performed  in  order  to  confirm  that  the 
conductivity  of  ice  is  too  low  for  the  ice-ice  mechanism  to  work, 
and  to  show  that  even  one  in  a  thousand  droplets  splashing  off  a 
riming  target  is  insufficient  to  explain  thunderstorm  electrifi¬ 
cation. 

Hailstones  which  occur  in  thunderstorms  carry  charges  up  to 
-100  pC.  The  simulated  hailstones  used  here  and  elsewhere  have 
all  been  conditionally  neutral  because  the  charge  acquired  by  the 
target  flows  to  earth  via  the  amplifier  and  chart  recorder.  It 
would  be  useful  to  place  a  charge  on  the  target  of  about  -100  pC 
and  then  rotate  the  target  through  a  mixed  cloud  to  determine  the 
effect  of  the  hailstone  charge  on  the  charging  process;  it  has 
been  suggested  that  this  may  limit  the  charge  transferred. 

The  work  involving  large  ice  crystals  is  continuing  at 
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APPENDIX  I 


THE  REPLICATION  OF  CLOUD  PARTICLES 

A. 1  The  Continuous  Pormvar  Replicator 

The  use  of  the  replicator  has  been  described  in  detail  by 
Hallett  (1976a, b)  and  Jayaratne  (1981).  The  device  is  used  for 
making  replicas  of  the  cloud  particles,  be  they  water  droplets  or 
ice  crystals.  The  cloud  particles  are  drawn  past  a  moving  film 
which  is  coated  with  formvar  solution  and  they  are  trapped  by  the 
sticky  surface.  The  formvar  dries  and  the  particles  evaporate 
leaving  a  plastic  replica. 

The  replicator  consists  of  a  film  support  arm,  take-up  reel 
and  the  drive  motor.  The  support  arm  (of  length  60  an)  fits  into 
a  smooth  metal  tube  of  the  same  length,  which  then  protrudes  into 
the  cloud  chamber.  The  transparent  16mm  film  was  transported  by 
two  dc  motors;  one  turns  the  take-up  reel,  and  the  other  is 
attached  to  five  of  the  rollers  mounted  on  top  of  the  replicator. 
The  formvar  solution  was  pumped  through  a  metal  tube  to  the 
applicator  head  at  the  end  of  the  support  arm.  As  the  film  moves 
past  the  applicator  it  is  covered  in  a  fine  coat  of  formvar,  the 
thickness  of  which  could  be  controlled. 

At  the  side  of  the  outer  metal  tube,  there  is  a  rectangular 
slot  of  dimensions  5.5  x  2.0mm.  The  end  of  the  tube  is  connected 
to  a  suction  pump.  When  the  support  arm  is  in  the  tube  the  slot 
is  located  just  after  the  applicator  head.  The  cloud  is  drawn 
past  the  film  and  impacts  on  it.  The  film  must  be  dry  before 
being  taken  up  by  the  reel,  and  so  the  film  is  moved  through  the 
rollers  on  top  of  the  replicator.  In  order  to  speed  up  this 
process,  cold  dry  air  is  blown  across  the  film.  The  flow  speed 
was  measured  with  a  flowmeter  and  was  normally  4  m  s  The  slot 
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was  heated  with  a  current  element  to  prevent  rime  building  up  and 
obstructing  tne  flow. 

Water  resistant  pens  were  used  to  mark  the  film  at  certain 
points.  Thus  the  information  on  the  film  could  be  correlated 
with  the  charging  current  shown  on  the  chart  recorders. 
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APPENDIX  II 


MEASUREMENTS  OF  THE  SURFACE  TEMPERATURE  OF 
A  RIMING  HAILSTONE 


Introd action 

A  platinum  resistance  temperature  sensor  was  used  to  measure 
the  surface  temperature  of  a  riming  cylindrical  target.  The 
sensor  is  the  same  diameter  as  the  simulated  hailstone  (5  mm)  and 
will  experience  the  same  effective  liquid  water  content;  thus 
the  need  to  calculate  the  mean  collection  efficiency  is  removed. 
Obviously,  if  the  sensor  is  thinner  than  the  target  it  will 
collect  more  of  the  droplets  due  to  its  higher  collection  effi¬ 
ciency,  and  so  would  over-estimate  the  hailstone  surface  tempera¬ 
ture.  One  of  the  first  attempts  to  predict  the  surface  tempera¬ 
ture  under  certain  conditions  was  made  by  Ludlam  (1951),  and 
later  extended  by  Macklin  and  Payne  (1967).  In  both  formulations 
the  surface  temperature  is  determined  by  a  balance  between  the 
rate  at  which  heat  is  released  on  freezing  and  the  rate  at  which 
the  heat  can  be  liberated  to  the  environment  by  forced  convection 
and  evaporation. 


The  Ludlam  Equation 


Ludlam  (1951)  formulated  an  equation  to  calculate  the  sur¬ 
face  temperature  of  a  riming  hailstone  as  it  moves  through  a 
supercooled  cloud.  The  equation  deals  with  the  rate  of  heat  loss 


by  convection  and  evaporation  and  heat  transfer  to  the  environ¬ 
ment.  The  rime  surface  temperature,  RT,  is  given  by: 

RT  =  AJKCT  +Lkdp)  +  B{CT  -  Lf)/(Ak  +  Bc^) 


where  A  =  ( 0 . 24 ( Re ) 0  *  6 )/d,  B 


Emv/V,  cw  =  specific  heat  of  ice, 
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CT  =  cloud  temperature,  d  =  cylinder  diameter,  dp  =  difference  in 
vapor  densities  at  the  cylinder  surface  and  the  environment,  E  = 
collection  efficiency,  K  =  thermal  conductivity,  k  =  coefficient 
of  diffusion  of  water  vapor,  h  =  latent  heat  of  evaporation, 

=  latent  heat  of  fusion,  m  =  liquid  water  content,  v  =  impact 
velocity . 


The  Macklin  and  Payne  Equation 

Macklin  and  Payne  (1967)  also  calculated  the  rime  surface 
temperature  using  a  similar  method  and  included  terms  to  account 
for  the  forced  convection  and  evaporation.  Their  equation  takes 
the  form: 


A{Pr1/3kTa  - 
B{Lf  +  Cw(Ta 
APr1/3k  +  Bci 


sc1/3lvd<ps 


Pe>> 


-  T  )  +  c.T } 
m  1  m 


Ts  =  ( Xf  +  X2)/X3 

where  A  =  (  0 . 24  (Re  ) 0 ' 6  )/d,  B  =  EWv/-j,  ,  =  specific  heat  of  ice, 

cw  =  specific  heat  of  water,  d  =  cylinder  diameter,  D  = 
coefficient  of  diffusion  of  water  vapor,  E  =  collection 
efficiency,  k  =  thermal  conductivity  of  air,  Lj  =  latent  heat  of 
fusion  of  water,  Lv  =  latent  heat  of  vaporisation  of  ice,  Pr  = 
Prandtl  number:  Pr  =  V  /K  ,  po  =  density  of  water  vapor  in 
environment,  ps  =  density  of  water  vapor  at  surface.  Sc  = 
Schmidt  number:  Sc  =  V  /D,  T  =  ambient  temperature,  T  =  time 
surface  temperature,  Tm  melting  temperature,  v  =  airspeed,  W  = 
liquid  water  content,  y  =  kinematic  viscosity  of  air. 
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APPENDIX  III 


THE  GROWTH  RATE  OF  CLOUD  HYDRGMETEORS 


Mason  (1971)  states  that  the  general  equation  for  the  rate 
of  increase  in  ice  crystal  mass  by  vapor  diffusion  takes  the 
form: 

dm  =  4TT  CD(  P  , ,  -  P  ,  )  (  A.  1  ) 

-T—  I  W  I  1 


where  the  constant  C  depends  on  the  shape  of  the  crystal,  D  is 
the  diffusion  coefficient  of  water  vapor  in  air  and  ^  and  p 

are  the  vapor  densities  at  the  crystal  surface  and  the  cloud 
respectively . 

For  a  hexagonal  plate,  C  has  the  value  of  2r/ir  if  it  is 

considered  to  be  a  circular  disc.  At  -  10"C,  ~  >  is 

0.21  g  m-3  and  the  rate  of  increase  in  mass  of  a  30  pm  diameter 

—  8  -1 

plate  is  approximately  10  g  s  The  rate  of  increase  of  mass 

per  unit  area  may  be  calculated  from  the  original  area  if  dr«r 
and  is  approximately  4  x  10-3  g  am-3  s-2. 

There  are  two  contributions  to  the  growth  rate  of  the  hail¬ 
stone:  firstly,  the  vapor  flux  to  the  surface  from  the  deep 
field  vapor  supply  and,  secondly,  the  vapor  flux  from  a  single 
droplet  freezing  or.  the  hailstone  surface.  These  two  contribu¬ 
tions  can  be  calculated  from  the  equations  given  by  Baker  et  al. 


(1987  )  . 

(i)  The  deep  field  vapor  flux  to  the  surface  is  given  by: 
vf  =  fD  (  p  w-  P  ^ /R  g  m-2  s'1 

where  f  is  the  ventilation  coefficient.  Thus  if  f=15,  d=2xl0  3 
m2  s  3,  R=2.5xl0  3  m  and  <p  w~  ^)=0.2  g  m  3  then  =2.4x10  3  mg 
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(li)  The  mass  growth  rate  due  to  tne  vapor  arriving  at  the 

surface  from  a  single  freezing  droplet  is  given  by: 

dm  =  3  TT  fDa  (p  -  f  .  )  g  s_i 
dt  °  1 

The  mass  to  droplet  lost  as  vapor  is  all  assumed  to  deposit  on 

the  hailstone  surface  and  is  given  by: 

M  =  3TT  fDa(  P  -  P  . )  t,  g 
1  o  1  1  t 

where  t^  is  the  freezing  time. 

Typically:  f=15,  d=2xl0"5  m2  s"1,  (  P  -  p 4  g  m~3, 
a=5  urn  tj=l  ms  and  so  M=5.6xl0  g. 

The  target  used  by  Jayaratne  et  al.  was  rotated  at  3  m  s  3 
in  a  cloud  of  droplets  of  diameter  5  pm  and  concentration  1000 
cc~3.  From  the  present  work  the  collection  efficiency  of  the 
target  for  this  droplet  diameter  is  0.2%  and  the  number  of 

droplets  swept  out  per  unit  area  per  second  is  300x1000x0.002  = 

-2  -1 
600  cm  s 

Thus  the  available  mass  fluxjt^f  ^s: 

=  5 . 6xl0-11  x  6x10 2  =  3 . 36xl0-8  g  cm-2  s_1 


The  total  rate  of  increase  of  mass  per  unit  area  of  the 
hailstone  is  just  the  sum  of  ana  ^2*  Thus, 

dm  =  2.43x10  ^  g  an  2  s  ^ 
dt 

The  ratio  of  the  rate  of  increase  of  crystal  mass  to  that  of 
the  hailstone  is: 


4x10 


2 . 43x10 


16.44 


that  is,  the  ice  crystals  are  growing  some  sixteen  times  faster 
than  the  hailstone  in  this  situation.  In  any  case  the  collection 
of  droplets  in  the  range5-8  pm  range  due  to  the  rime  feathers 


64 


still  does  not  provide  as  much  vapor  as  from  the  droplets  in  the 


cloud.  Therefore,  the  conclusion  of  Jayaratne  et  al.,  tnat  the 
charge  transfer  in  these  experiments  is  controlled  by  the  vapor 
diffusion  to  the  target  surface  from  droplets  in  the  cloud  still 
stands  . 
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APPENDIX  IV 


THE  EVENT  PROBABILITY  OF  A  CYLINDRICAL  TARGET 


The  collection  efficiency  of  a  cylindrical  target  for  super¬ 
cooled  water  droplets,  can  be  considered  to  be  the  same  as  the 
collision  efficiency  as  only  a  few  droplets,  1  in  1000, 
Aufdermaur  and  Johnson  (1972),  will  rebound;  the  rest  will 
freeze  on  impact.  For  ice  crystals,  however,  the  collection 
efficiency  is  not  the  same  as  the  collision  efficiency,  as  the 
ice  crystals  have  a  finite  separation  probability  which  increases 
with  crystal  size  and  impact  velocity.  The  event  probability  is 
defined  as  the  probability  of  a  crystal  impacting  on  a  target  and 
rebounding.  Therefore,  the  event  probability  is  a  product  of  the 
collision  efficiency  and  the  separation  probability.  The  colli¬ 
sion  efficiency  equations  of  Ranz  and  Wong  are  not  applicable  to 
ice  crystals  due  to  their  different  shapes  and  densities. 
Jayaratne  (1981)  states  that  the  event  probability  of  ice  crys¬ 
tals  is  directly  proportional  to  the  total  charge  separated  in  a 
given  time. 

The  simulated  hailstone  targets  consist  of  a  rod  and  wire  of 

diameters  R  (5.0  mm)  and  D„  (0.5  mm)  respectively.  It  is 
r  w 

assumed  that  the  two  targets  experience  the  same  cloud  conditions 
at  all  times.  The  charging  currents  to  the  roc  and  the  wire  are 
I_  and  L.  respectively.  The  volume  of  cloud  swept  out  by  either 
target  is  directly  proportional  to  the  target  diameter  and  the 
charging  current  is  directly  proportional  to  the  event 
probability.  If  the  event  probabilities  of  the  two  targets  are 
Er  and  respectively,  then: 


and  therefore: 


erdr/ewdw  =  V1* 
**  a  'Wr/Vw 
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The  wire  diameter  is  small  compared  to  the  rod  diameter,  and  so 


it  is  assumed  that  it  has  an  event  probability  of  unity,  and  if 
the  diameters  are  substituted,  then: 

E*  =  °-5  V5  *W 

This  relationship  is  valid  for  all  crystal  sizes  and  impact 
velocities  (except  very  low  velocities),  because  is  assumed  to 
be  unity. 

The  values  of  the  event  probability  given  by  the  ratio  of 
the  charging  currents  were  confirmed  by  the  experiments  described 
in  Section  5. 
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